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ABSTRACT 


An instrument to measure atmospheric pressure at the Earth's surface 
from an orbiting satellite would be a valuable addition to the expanding 
Inventory of remote sensors. The subject of this report is to describe 
such an Instrument - the Microwave Pressure Sounder (MPS) - based on an 
active rillimeter wave technique. It is shown that global-ocean coverage 
is attainable with sufficient accuracy, resolution and observational fre- 
quency for meteorological, oceanographic and climate research applications. 

Surface pressure can be deduced from a measurement of the absorption 
by an atmospheric column at a frequency in the wing of the oxygen band 
centered on 60 GHz. An active multifrequency instrument is needed to make 
this measurement with sufficient accuracy. The selection of optimum 
operating frequencies is based upon accepted models of surface reflection, 
oxygen, water vapor and cloud absorption. By combining the ratios of 
transmission of three pairs of frequencies the measurement is made es- 
sentially independent of variations in the atmospheric temperature and 
water vapor profiles, cloud cover and sea state. 

Numerical simulations using a range of real atmospheres defined by 
radiosonde observations have been used to validate the frequency selection 
procedure. This also shows that the error introduced by unmodellable 
atmospheric variations is less than 0.4mb. Other sources of error are the 
instrumental noise and the statistical nature of the ocean surface re- 
flectivity. Analyses are presented of alternative system configurations 
that define the balance between accuracy and achievable resolution. An 
MPS system which requires minimum development of new components, when 
used in an 800 km orbit, can provide measurements of surface pressure 
with a predicted r.m.s. error of 1.2 mb and a surface resolution of 
10 X 100 km. 

■ 

The global coverage of pressure measurement provided by the MPS has 
been examined for a variety of implementation schemes. Oceanic coverage 
which is comparable to the existing density of data points on the continental 
land masses can be achieved with a network of three instruments in sun 
synchronous, near-polar orbits. 

An ongoing U.K. - U.S. development program is being directed 
towards an early proof-of-concept flight for the MPS on a Shuttle/Spacelab 
sortie mission. Subsequent operational deployment of the MPS will provide 
valuable data for numerical weather forecasting and oceanography while the 
accumulation over a long time period of consistent measurements of the global 
pressure field will contribute significantly to the data resources needed for 
climate studies. 


- 2 - 




1 , M etoorolof^lcal Backg round 

Global meaRuromenta of the Earth’s atmospheric prcsnuro are important 
In synoptic meteorolop,y , numerical weather forecasting, atmospheric dynamics 
and climate studies. At the present time pressure data are gathered prin- 
cipally from land based weather monitoring stations, and arc supplemented 
over the oceans by reports from ships and aircraft. Stations reporting 
surface pressure to the National Meteorological Centre at OOZ on February 1, 
1976 are shown in Fig. 1 and it is apparent that over the tropical and 
southern ocean areas in particular there is a deficiency of data. Such non- 
uniform coverage is inadequate in terms of the requirements for a set of 
global meteorological data sufficiently accurate to provide the initial 
conditions for a numerical forecast of the weather. 

The World Meteorological Organization has specified (WMO, 1973) a set 
of observational requirements for the First GARP Global Experiment of the 
Global Atmospheric Research Programme (CARP) . Measurements of pressure in 
data-sparse regions are required with a horizontal resolution of 500 km and 
an accuracy of ± 0.3% (equivalent to ± 3 mb at the surface) . However, a 
recent survey of user needs (JPL, 1976) indicates that an accuracy of 
± (1-2) mb in sea level pressure is desirable and it is probable that a 
± 3 mb specification is near the limit of usefulness. Resolution require- 
ments ranged from 250 km down to 10 m but the lower figure is unrealistic 
and not consistent with the accuracy requirement. 

The cost of acquiring this quantity of data by conventional land and 
ocean monitoring stations would be unreasonable (Mason, 1970) and the 
variation of instrumentation and calibration between stations would reduce 
the reliability of the measurements. In contrast, the provision of 
operational satellite borne instruments for remotely sensing atmospheric 
temperature has proved inexpensive and has given accurate consistent 
global coverage for several years. Remote sensing of atmospheric water 
vapor has been demonstrated with the Nimbus E Microwave Spectrometer 
(Staelin et al., 1976) and Seasat A, due to be launched In 1978 will carry 
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an instrument to obtain surface winds from sea state measurements. The 
development of an Instrument to sense pressure from a satellite ol at form 
would be another major advance towards the WMO requirements. It would seem 
appropriate to aim for an Instrument with an accuracy of .i(l-2) mb and 
capable of a horis'-ontal prld spaclnp, of 100-2()f km to i < respond witli the 
data at present available fr(*m existinp, land stations. 

A network of geostationary natellJtes at !lb,000 km altitude provides 
near global coverage with a minimum number of satellites. Ilovever , sat<’l- 
lites in near-polar orbits at 1000 km .lUitiide. arc required to obtain data 
for high latitudes and to allow microwave sensing with su' C'cicnt r'*sol ition. 
Both types of satellites will be elements in any operational systei. .or com- 
prehensive monitoring of meteorological parameters within the fo v liable 
future. The active microwave technique for pressure i sur l o be dis- 
cussed here is only feasible from the love^* .fered by a polar orl iter. 

The accuracy of 10.3% in pressure measurement corresponds to a height 
change of only 25 m at sea level so that the height at which the pressure 
is mea.sured must be known to within this value. This constraint effectively 
precludes at present the possibility of measuring the three dimensional 
pressure field. However, a knowledge of the surface pressure together with 
a measured temperature profile allows the hydrostatic approximation to be 
used to obtain the variation of pressure with height. An instrument to 
measure surface pressure is therefore entirely adequate when used in con- 
junction with a temperature sounder to define the mass field. 

The proposed method for measuring surface pressure has been specifi- 
cally designed for use over the oceans. However, measurements over land 
may give a valid indication of pressure, but only if a suitable average 
height can be used and if the surface has appropriate reflection character- 
istics. This is not a serious disadvantage since oceans cover some 70% of 
the Earth's surface and while the existing network of weather monitoring 
stations gives satisfactory coverage of most of the continental land masses, 
the oceans are almost totally neglected. The satellite surface pressure 
data for the oceans will complement present land data to approach the 
GARP requirement. 


Methods for the Remote Measurement of Pressure 


Atmospheric pressure may in principle be determined remotely from 
measurements of the characteristics of the propagation of electromagnetic 
radiation through the atmosphere. The accuracy requirement necessitates 
the use of an active system and several alternative methods have been 
proposed. The feasibility of their Instrumental implementations has been 
assessed in the literature (Singer 1968, Mason 1970, Smith et al 1972, 

AMWR 1975) and the following brief summary of some of the limitations of 
these methods indicates the reasons for adopting the technique to be des- 
cribed in this report. 

The three dimensional pressure field can in theory be determined by 
a range— gated radar or lidar instrument. The pressure dependent contri- 
butions to atmospheric backscatter are too weak to be measured but a dif- 
ferential absorption, 0.76 ym lidar relying on aerosol or molecular Rayleigh 
backscatter appears to be feasible with probable technological advances 
within the next decade. However, to achieve the coverage and accuracy re- 
quirements it would need more primary power than is available in currently 
planned operational meteorological satellites. Also, because of the in- 
ability of infrared radiation to penetrate cloud its operation would be 
restricted to measurements down to the surface in clear skies and above 
cloud measurements elsewhere. 

Although the impracticality of pressure profile measurements leaves 
only methods that measure pressure at a reference level this is not a 
serious limitation. Such data in conjunction with measured temperature 
the full pressure field to be deduced using the assumption 
of hydrostatic equilibrium. For this purpose, the surface is the most 
satisfactory reference level. The limitations on coverage in cloudy con- 
ditions with infrared instrumentation also apply to visible techniques 
but to a much lesser extent at microwave frequencies. It is apparent that 
mtCAowttve techniques for measuring AuA^ace pressure offer the best possi- 
bilities. 

Transmissivity measurements through the limb of the atmosphere can 
provide pressure at a chosen reference level. Above 12 km the pressure 


can be obtained by passive nitr.rownve measurement^ but Inadequate knov/lcdp.c of 
the temperature and water vapor profiles would not allow surface pressure to 
be deduced from this to sufficient accuracy. To obtain pressure at a lower 
altitude from a limb measurement would require an active occultation techni- 
que with the inconvenience of two satellites whose positions must be con- 
tinuously monitored to a high accuracy. The horizontal resolution with this 
method is several hundred kilometers and is thus insufficient to resolve some 
of the most interesting features of the pressure field. 

Similar limitations affect bistatic refractivity measurements re- 
cently proposed to determine surface pressure by reflecting a microwave 
signal from the ocean surface at a low grazing angle. In this case accurate 
knowledge is required of the water vapor field and the geo id. The 3 mb 
accuracy and 500 km horizontal separation of measurements can only be ap- 
proached with a constellation of high-altitude inter-communicating satellites, 
the orbital characteristics of which are incompatible with other demands on 
a meteorological satellite network. 

An optical path length method measuring the dispersion between two 
different microwave frequencies requires a time resolution much less than 
the length of a pulse reflected from the rough ocean surface and is there- 
fore not practicable. 

Differential absorption measurements in and near the millimeter wave- 
length absorption band of oxygen by reflecting signals from the Earth's sur- 
face provide the basis of the chosen method. It will be shown that measure- 
ments at several frequencies in the wings of the oxygen band centered on 
60 GHz combined with measurements closer to the water vapor line at 22 GHz 
have the potential of achieving the required 1-2 mb accuracy of surface 
pressure determination without demanding an accurate knowledge of the tempera- 
ture and water vapor fields. Spatial resolution of the order of 200 km is 
possible and coverage while not completely global, is better than for the 
other methods considered, furthermore the instrument is suitable for flight 
on operational meteorological satellites and does not demand excessive pri- 
mary power. 


3. An Outline of fhe Microwave Pressure Sounder Concept 


In this section the conceptual basis of the proposed technique is out- 
lined and its plausibility established. An introduction is presented to the 
various physical aspects of the problem which together will be used as the 
basis for the Instrument design criteria. The necessary details of the 
analysis will be given in subsequent sections but the symbols to be used are 
defined here and an indication is given of the order of magnitude of the 
principal quantities. 

3.1 Principle of the Method 

With a Microwave Pressure Sounder (MPS) in orbit as depicted in Figure 
2 the surface pressure is deduced from measurements of the ratios of atmo- 
spheric transmissivity at pairs of frequencies. The primary pair of pressure 
sensing frequencies lie on the lower wing of the broad 60 GHz oxygen absorption 
band, one sufficiently far into the band to be significantly affected by changes 
in the total mass of oxygen in the transmission path with the other frequency 
on the extreme wing of the band and therefore affected very little by oxygen 
absorption. The ratio of transmissivities, or differential absorption, is a 
measure of the total oxygen in the path and hence on the partial pressure of 
oxygen at the surface. Since oxygen is a uniformly mixed constituent of dry 
air this is also a measure of the total surface pressure. Differential ab- 
sorption measurements are necessary at two other pairs of frequencies to 
make the pressure measurement insensitive to the atmospheric temperature pro- 
file, water vapor content, cloud cover and sea state. To make these absorption 
measurements an active system is required. For each pair of frequencies, 
signals of equalized power are transmitted and, following reflection from the 
ocean surface, the ratio of the returned echo powers is measured. When the 
operating frequencies have been chosen correctly, a combination of the three 
measured ratios is dependent only on the surface pressure. Thus, the output 
of the MPS is simply related to the surface pressure and does not rely on sub- 
sidiary meteorological measurements for interpretation. 

3.2 The Reflected Signal 

The radar equation generally used to give the magnitude of the signal re- 
ceived at the antenna from an extended reflector such as the sea surface is 
(Skolnik 1970). 






3.1 


Pr - 


2 2 

Pj Saif's^ T (0) a (6) 
(4tt)^ R^(6) 


where the symbols have the following meanings and are identified on Fig. 3. 


Pr received power at antenna 

Pfj, transmitted power 

G antenna gain 

03<l>3 antenna 3db angular beam widths 

X signal wavelength 

a(0) radar cross-section of reflector 

t( 6) total atmospheric transmission coefficient (one-way) 

R(0) distance between receiver and reflector 

6 angle from nadir 


The 3db beam widths are related to the antenna gain thus 

03(j)3 s 4 ^/q ^ ^ 

and the gain is related to the area. A, of the antenna thus 

G « 4ttA/x2 2 3 

It is convenient and less confusing to restrict the equation to the nadir 
case (0«O) which is of principal concern. It can be readily extended as 
necessary for the near nadir situation which will be considered later. 

Using the relationships 3.2 and 3.3 the radar equation may then be written 
in the simpler form 


where h is the satellite altitude. 


3.4 


Transmitter 
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r Atmospheric 
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o Backscatter cross-section 


Figure 3, Symbol Definition 
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The significance of the various factors in equation 3. A can be readily 

seen 

: the power reaching the sea surface is the transmitted power 

reduced by the atmospheric transmissivity factor 

; this is a geometrical factor giving the fraction received by 
an area A at a range H from an isotropic emitter. 

o ; represents the reflecting ability of the sea surface in re- 

lation to an isotropic reflector and includes the directionality 
of the reflection as well as the reflectivity of t)ie surface. 

T ; a further reduction due to the return path through the atmosphere. 

For the purpose of instrument design it is better to put the equation 
in terms of the radius of the antenna, r, using A - Trr^ and to Introduce 
the efficiencies, n.j,, of the transmit and receive antennae feed systems. 

In addition the noise, N, must be specifically included since this has a 
dominating influence on the system design. Introducing the detected power 
Pjj. which is the sum of the detected return echo. P^. and the noise we may ’write 
Pd = Pr + N - 0.25 ^ ,2 ^0 + N 3.5 

The potential of the method for achieving the required accuracy can be 
demonstrated by estimating the signal to noise ratio using sensible values 
this equa-ion. For a IW source, a .5 m diameter antenna with efficiencies 
n.j,. = 0.5 at a height of 1000 km above a sea surface of radar cross-section 

10 db and an operating frequency such that t - 0.5 the received signal is 
'V'lO” W. The noise equivalent power is 


N » 4kT (Af/t) F - , 

Ja O 

where k is Boltzman’s constant and 

T receiver temperature 300 K 

Af the IF bandwidth lOOkHz 

t integration time 10 s 

F receiver noise factor 10 db 


These values give a S:N ratio of 600:1 which Is compatible with the 0.3Z 
accuracy requirement. 

3.3 Design Parameters 

The physical factors influencing the magnitude of the received signal 
can be clearly seen from the various terms In the radar equation (3.5). In 
addition the parameters avallabe for system design can now be Identified. 

The total atmospheric transmissivity should first be separated Into terms 
covering oxygen, water vapor, t^, and other atmospheric constituents, 

'a- 

Pd * 0.25P^n.^nR ^ t2 t2 ,2^0 + n 3.7 

P^ and Pjj are the source and detected powers. Source power of the order 
of IW is required and solid state Impatt oscillators are probably the most 
suitable for space applications because of their low voltage requirements 
and high efficiency. The power available from single devices has an upper 
limit of about IW but by 1978 sources delivering several watts should be 
available when techniques for combining several diodes together are ex- 
tended to frequencies up to 70 GHz. 

The antenna feed efficiency, n, has a value dependent on antenna 
design and system losses, but a working figure of 0.5 can be used. The 
antenna radius, can be chosen to optimize system performance by balancing 
slgnal-to-nolse error, which demands a larger antenna, against a statis- 
tical error due to sea surface effects, which demands a smaller antenna. 

The optimum value for r is of the order of 25 cm so that antenna subsystem 
design problems due to excessive size are not encountered. The field of 
view at the surface is determined by r and it has a diameter of 30 km 
for r ■* 0.25 m with the satellite at a height of 1000 km. 

The orbit height h should preferably remain reasonably constant to 
avoid the need for a receiver with a large dynamic range. A high orbit 
provides more favorable geometry for meaningful measurements away from 
the sub satellite path but as h is decreased the signal-to-noise ratio 
and consequently the accuracy of the measurement are improved. Too low 
an orbit however will limit the operational lifetime of the satellite and 
a value for h between 250 km and 1000 km is desirable. 


*, 
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Tq the transmissivity of the atmosphere due to oxygen over a double 
path can be written as 10 "° *^^0 where Tq is the integrated absorption 
coefficient in db for a single path. Fq can be calculated from theoretical 
models of the oxygen absorption band (Waters 1976, Rosenkranz, 1975) in 
terms of the variables surface pressure, pg, and temperature profile, T(z) , 
for any value of the disposable parameter, v, frequency. The zenith 
attenuation due to atmospheric oxygen as a function of frequency is shown 
for a U.S. standard atmosphere in Fig. 4. The principal operating fre- 
quency of the instrument should be chosen around 52-53 GHz. Further into 
the band the signal and S:N ratio are reduced while at lower frequencies 
the effect of surface pressure on the signal is reduced. 

The presence of water vapor in the atmosphere has both indirect and 
direct effects. It affects Tq because it modifies the pressure broadening 
of the lines within the oxygen band. It contributes to absorption in the 
spectral band of interest through which can be written as 10"°*^^w 
and Fy, the integrated absorption coefficient for water vapor, is cal- 
culable from existing models in terms of the pressure profile of water 
vapor ^(z) and T(z) for given values of v. Finally, water vapor can 
make a significant contribution to the total surface atmospheric pressure. 
The water vapor spectrum is also shown in Fig. 4 for a surface density 
of 12.3 gm m“^, decaying exponentially with a scale height of 2,2 km. 

With the exception of ozone, other atmospheric gases do not signi- 
ficantly absorb in this region of the spectrum, but clouds and rain both 
absorb and scatter. Numerous weak ozone lines exist throughout the 
microwave spectrum but only close to the line centers is the absorption 
strong enough to affect this experiment. These frequencies can be easily 
avoided . 

The frequency dependence of atmospheric constituents other than 
oxygen and water vapor can be modelled by an absorption coefficient, 
which is a second order function of frequency, viz 

F^ ■ a + bv + cv^ 3.8 
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The coefficients a, b and c are highly variable both in space and time 
because equation 3.8 represents the effects of cloud and rain. The 
experiment must be designed so that the transmissivity measurements allow 
these coefficients to be cither evaluated or eliminated when deducing 
surface pressure. 


The radar cross-section of the sea surface, o(0), depends upon the 
sea state, the angle of incidence and the frequency. The smooth variation 
with frequency is adequately modelled by a second order function so that 
equation 3.8 also Includes sea surface frequency dependence. o(0) has a 
maximum value, o°, for normal incidence. It reduces more rapidly for 
calm seas than for rough seas as the angle of incidence decreases. For 
most sea states it should be possible to make measurements with the pro- 
posed instrument at angles of up to about 15° from the vertical. Such 
measurements away from the subsatellite path will allow the surface 
pressure gradient to be determined and this is of much greater value than 
pressure only along the orbit path. For a 1000 km orbit altitude the 
15° limitation means that the swath width at the surface is about 500 km 
across which significant pressure differences frequently occur. Variations 
in the sea state due to surface winds can change the average value of 
0° by a factor of 10 or more so the dynamic range of the instrument must 
accommodate this as well as atmospheric absorption variability. 


Perhaps the most important effect of the sea surface is that because 
of its constantly changing nature, both spatially and in time, the co- 
herence of the reflected signal Is destroyed. The Intensity of the return 
echo varies statistically with a Rayleigh distribution, the standard 
deviation of which is equal to the mean. The average of a large number 
of independent samples is required to achieve the necessary accuracy. 

The relationship between the sample collection rate and the antenna radius 

is discussed in detail in Section 11 together with the consequences this 
has for system design. 

The noise N Includes contributions from background emission, atmospheric 
reflection and instrumental noise. It depends most importantly on the 
bandwidth, Af . of the detection system and on the integration time. t. 
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Cholce of t will, for a given satellite speed, determine the ground resolution 
of the instrument along the orbital path. A typicnl satellite speed is 
'^7kras~^ so that the 500 km resolution required for GARP can be achieved 
with t up to 1 min. However, shorter integration times and consequently 
better spatial resolution are anticipated because although 500 km is ade- 
quate resolution for the input data set of numerical global circulation 
models many meteorologically significant features in the pressure field have 
dimensions leas than this. 

Backscatter from atmospheric constituents is small except perhaps 
from heavy rain, but under these conditions the increased absorption will 
probably make the instrument unusable. To reduce noise to a minimum the 
bandwidth should be limited to that necessary to include all the return 
signal. A direct measurement of noise must be made to enable the best 
estimate of signal level to be determined. 

3.4 Summary 

This overview of the proposed technique for determining atmospheric 
pressure at the surface of the Earth from transmissivity measurements has 
Identified the following parameters, the values of which are at our dis- 
posal for designing a system with optimized performance: 

transmitter power 
h orbit height 

r antenna radius 

V operating frequencies 

t integration time 

4f receiver bandwidth 

Measurements must be made at several frequencies to enable signal 
level changes due to pressure variation to be distinguished from the 
effects of all of the other factors mentioned here. It is necessary 
next to discuss in detail the theory of these physical factors so as to 
establish a basis for system design and choice of operating frequencies 




uaing the disposable parameters that have been identified. It will bo 
shown that an instrument can be designed for which 

• the sensitivity of the measurement of pressure to the 
temperature profile T(z) ard to the presence of water 

vapor in the atmosphere is reduced to an insignificant 
level; 

• multifrequency measurements are used to eliminate the 
effects of the other unknown quantities which are 
described by the variable coefficients a, b and c; 

• averaging over a large number of independent samples 
provides the accuracy demanded ; 

• useful near global-ocean coverage can be obtained; 

• available technology is used for all components except 
the transmitter source and these are currently being 
developed. 
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A . Absorption of Mlillmotor Wavon by Oxyg on 

To be able to determine pressure from multifrequency me.nsurementR of zenith 
attenuations detailed knowledge must be. available of the oxyp.en absorption co- 
efficient as a function of frequency, presnure and temperature. This nertlon 
auiranarl7.es the relevant characteristics of the most recent theoretical apectro- 
acoplc model of oxygen (Rosenkranx 1975). The theory produces results In close 
agreement with experimental data (Llebe et al. 1977). 


A.l Absorption Coefficient 

The transmission of radiation at a frequency v through an ab.sorbing 
msdlum is generally described by the equation 

Iy(x) » ly(0) exp (-Y(v)x) A.l 

where the Intensity Iv(x) after the radiation has travelled a distance x is 
given in terms of the initial intensity lv(0) and the absorption coefficient 
y(v). 


For a path through an inhomogeneous medium the absorption can be obtained 
by integrating y('o) . In particular for a vertical path through the atmesphere 
the absorption of millimeter waves due to oxygen is expressed in terms of the 
integrated absorption coefficient To given by 

To ^*2 

where z is the vertical coordinate. 

A. 2 The Microwave Spectrum of Oxygen 

The absorption by molecular oxygen of electro-magnetic waves in the 
millimeter spectral region, was first described quantitatively by Van Vleck 
(19A7). The quantum number N, of which only odd values are allowed, describes 
the rotational angular momentum of the molecule. This couples with electronic 
spin to form a triplet of states with total angular momentum J = N~l, N. or 
N + 1. Selection rules permit magnetic dipole transitions with AJ = ±1, 0. 

The transitions A J • +1 give rise to a broad band of lines between 30 and 
70 GHz and an Isolated line at 118 GHz. The resonant frequencies are denoted 









: > ■ 


Vjj or for transitions (J » N) -*■ (J = NHl) and (J - N) -► (J - N-i) re- 
spectively. AJ » 0 transitions correspond to "non-resonant” absorption at 
zero frequency. The total absorption at any frequency is calculated by 
sutnming over those lines with significant absorption and it is found to be 
sufficient to Include only those lines with N < 39. 


In the pressure range from 1 to 1000 mb the dominant line broadening 
mechanism is intermolecular collisions. Rosenkranz gives the absorption in 
db km“^ as 


Y(v) » 1.434 p F(v) 


4.3 


where the pressure p is in mb, the frequency v is the GHz and the temperature 
T is in K. The function F(v) is a line shape factor which together with 
describes the spectral intensity distribution and is given by 


F(v) = p 




N 


N 




where 


f*(v) = vn) 7n 

(v - v^)2 + (pwjj)2 


4.4 


4.5 


It can be seen that the summation is over a set of lines at both positive 
and negative frequencies with line shapes (eqn. 4.5) of Lorentzian form 
modified by the interference coefficients y^. The final term of equation 4.4 

is the contribution from the zero frequency lines where w^ is the non-resonant 
linewldth. 


The population factor is given by 

2n + 1 


“ 0TT25T (-2.0685 N (N+1)T) 


4.6 


and the line an^lltudes dH are 

N 


N 


'N 


N(2N + 3) 
(N+1)(2N+1) 


\h 


(N-fl)(2N-l) 
N(2N + 1) 


4.7 


4.8 


Following Rosenkranz the interference coefficients are approximated by 


2d: 


N+2 


w, 


N 


^‘^N-2 ”n 


N ^ ± ± 

v" - V v^. ” ' 

* N N+2 


'N 


N-2 


_k 

± 

^N 


vi +60 

N 


4.9 


where wl and wl are near-diagonal elements of the [w] matrix when 2xp[w] is 
N N 

the transition rate matrix. The approximation restricts the coupling of each 

positive frequency line to the two adjacent lines in the same branchy to the 

line at zero frequency and to the lines at negative frequencies. Two con- 

t 4* 

ditions can be applied to the elements w^^ and w^^ which then allow them to be 
evaluated by an iterative procedure from the relationships 


w. 


N 


''b “ ”n “ 


4.10 


and 


N-2 


w, 


N 


*N^*N -2 


4.11 


f 

The procedure is initiated by arbitrarily setting = 0 . This produces a 
relatively large error in the Interference coefficient for large N but it is 
not significant since these are weak lines and the error is progressively 
damped by application of the conditions 4.10 and 4.11. 


The experimentally measured values of frequencies, and linewidths, 
w*, given in Table 1 were used to compute the microwave spectrum of oxygen. 
The non-resonant linewidth is the least accurately known parameter and it 
has a significant impact on the integrated absorption calculations away from 
the individual line centers and in the wings of the broad band. Rosenkranz 
uses the value 

w^ » 0.48 (-^) MHz mb" ^ 4.12 

which is the average of measurements made by Kaufman (1967) and Maryott and 
Birnbaum (1960). This value is preferred to a later measurement by Liebe 
et al (1977) which is 0.54 MHz mb”^ at 300 ”k because it gives a better fit 
to Reber*s (1972) opacity measurements from 66 to 71 GHz (Fig. 5) and biebe's 
technique is n. - particularly sensitive to variations in 



'■"t'lo i. Pr«,„™cu.s, I.<„e„ldtl,s of tl,.. Oxygen Klernwevo 

Sivjci I'lim (Kr* M.cba et al, 1977 ) 


Zenith Attenuation (db) 



Figure 5. Rosenkranz theory of oxygen absorption fitted to 
Reber's measurements of Zenith attenuation. 
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5. Microwave Absorption by Water Vapor 

As well as modifying the linewidths of the oxygen spectrum (Section 4.3) 
atmospheric water vapor affects the measurements to be made with the pressure 
sounder by absorbing millimeter waves and by contributing to the total 
surface pressure. 


The ifticrowave spectrum of water vapor is caused by electric dipole 
transitions between rotational states of the molecules. The lines of parti- 
cular interest here are centered at 22.235 and 183.310 GHz and the best 
agreement between theory and experiment is obtained by using a kinetic 
lineshape (Waters, 1976). However, in the window region calculated absorption 
is significantly less than that measured. Waters shows that for frequencies 
up to 100 GHz the discrepancy can be satisfactorily accounted for by grouping 
together the contributions from all lines other than that at 22 GHz and by 
including an empirical correction term. The absorption coefficient in db km"l 

Y^(v), is then 

/ . I 3.12 X 10^ exp(-644/T) ^ 012 I 5.1 

^w''* .j.3/2 j (22.235^ - v^) 2 + 4v^ Av^ J 


where the linewldth in GHz Is 

- - (ife) ) 


.0.626 


5.2 


The water vapor density p is in gm m“3 and all other symbols are as previously 
defined. Equations 5.1 and 5.2 enable the absorption due to atmospheric water 
vapor to be determined for given values of p, T, p, and v. These computations 
agree with measured absorption within the 5 to 10% accuracy of measurement 
generally achieved with this difficult test gas (e.g. Liebe and Welch, 1973). 
The anomalously high absorption away from the line centers has no satis- 
factory theoretical basis. At higher frequencies there are strong indications 
that water vapor dimer effects contribute significantly to the absorption 
(Harries and Ade 1972). However, below 100 GHz there is no evidence of a 
dimer associated spectral structure which could affect the choice of pressure 
sounder operating frequencies. 
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To account most conveniently for the contribution from water vapor to 
the total surface pressure, the effect of absorption due to any partial 
pressure of water vapor on the instrumental pressure index obtained from 
multifrequency measurements should be equivalent to the effect of an equal 
pressure of dry air. 
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6. Zenith Attenuation 
6.1 The Vertical Path 

The equations of Sections 4 and 5 permit the evaluation of the absorption 
of microwaves with frequencies up to 100 GHz by a homogeneous atmosphere of 
specified pressure, temperature and water vapor content. By integration, 
absorption through a vertical and thus non-homogeneous path can also be ob- 
tained. This has been computed for the U.S. Standard Atmosphere with added 
water vapor shown in Fig. 6. The water vapor content was specified in terms 
of its mass mixing ratio, S, such that at any pressure, p, 3 is given in 

terms of the surface mixing ratio 3<, and pressure p. by 

8 8 

6 - 6g 6.1 

Thus the water vapor falls exponentially with height with a scale height of 
about 2 km compared to the 8 km scale height for pressure. 

To carry out the numerical Integration the atmosphere must be divided 
into a reasonable number of layers. The absorption of each layer was approxi- 
mated by the value for an equally thick homogeneous layer using average 
values of temperature, pressure, and water vapor mixing ratio. Several dif- 
ferent methods have been used to evaluate these average values and various 
numbers of layers tried. It was found that as the number of layers was 
Increased beyond 24, divided equally between the troposphere and stratosphere, 
the resulting changes In calculated absorption were not significant for the 
purpose of system design. However, for consistency with the numerical simu- 
lations of Section 8 it Is convenient to use the heights and temperatures of 
the standard atmosphere at the 36 standard pressure levels between 1013 and 
2 mb used for reporting radiosonde observations. 

The zenith attenuation thus calculated Is shown in Figure 4 for a water 
vapor mass mixing ratio at the surface 3g ■ 0.01 which corresponds to a 
density of 12.9 gm m“^. 



Height (km) 
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6,2 Derivatives w.r.t.T, p and Water Vapor 

To enable the operating frequencies of the pressure sounder to be chosen 
80 that the measurement is insensitive to temperature profile and water 
vapor content of the atmosphere the change in integrated absorption coeffi- 
cient due to changes in T(z) and water vapor content must be determined. 

The sensitivity of the sounder to pressure can be obtained from the deriva- 
tive with respect to pressure. The method by which the derivatives are 
calculated is given here but the reasons for defining them in the following 
way are given in Section 7. 


For convenience we introduce a total integrated absorption coefficient • 
r.j,, which is the sum of the coefficients, for oxygen and for water 
vapor. At a particular frequency can be expressed as 


r + r 

0 w 


To 


\Po/ \'^of 


x(y+i) 


6.2 


where the surface pressure, Pg, and ten?>erature, Tg, are normalized to some 
standard values p^ and T^ for which the coefficient has the value The 

indices x and (y+1) are functions of frequency. Thus the derivatives w.r.t. 
fractional changes in the surface pressure and temperature are, respectively 


and 



6.3 


6.4 


The index x was determined by repeating the calculation of Section 6.1 
for an atmosphere with 2% increase in surface pressure and appropriate in- 
creases at other levels to maintain hydrostatic equilibrium. The values of 
zenith attenuation, and F^, thus give x 

X - 

log (1.02) 


The index (y+1) can be similarly obtained from a calculation of zenith 
attenuation for a 2X increase in temperature throughout the atmosphere so 
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that a change in T is taken to be representative of a change in the profile 
T(z). These are then used with equations 6.3 and 6. A to give the derivatives 
which are shown in Fig. 7 as functions of frequency. In Figures 8 and 9 
the data for the important wing regions is presented with an expanded fre- 
quency scale to show the fine structure associated with individual oxygen 
lines within the broad band. 

The effect of a change in the water vapor content was also determined 
by recalculating the zenith attentuation when the surface mixing ratio 6^ 
is increased by 2%. However it is plotted in Figure 7, 8 and 9 as the 
change, (6r^)^» in directly in db for effectively a 100% change in 
water vapor content calculated by multiplying the difference in for the 
2% change in 3g by 50. 


<1 
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Figure 7. Derivatives of the integrated absorption coefficient w.r.t. 
pressure, temperature and water vapor. 
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7. Choice of Operating Frequencies 


7.1 The Iixdex of Surface Pressure 


llie basic measurements from which surface pressure can be deduced are of 
the return echo signal strengths, P_, at several frequencies. We define an 
index of pressure, S, which is obtained by combining the measurements to- 
gether as a product thus 




and the indices r^^ are introduced to give greater flexibility in the choice 
of operating frequency. From equation 3.7 we can write the return signal at 
a frequency as 

^ (VtV')v^ ('o ("a 

SO that S can be expressed as the product of instrumental, atmospheric and 
background factors: 

S = S- S. S_ . 7.3 

I A B 


The instrumental factor 


, _ n 1 0.25 /_ , V )r^ 

>T “ n r— (P„ri„rir.r^ 1 i 

I i I h^ \ T T R / v^l 


should be a constant and this requirement affects the system design but not 
the choice of frequencies. Measurements will be made of the ratio of return 
signals at a pair of frequencies since experimentally this is easier than 
measuring absolute intensity. With this technique 

r^^ = -r^ ^ ^ (for odd values of i) 7.5 

and it can be seen that neither the antenna radius, if it is the same for all 
frequencies, nor the orbit height will affect Sj. If the transmitted 
powers at each pair of frequencies are controlled to have a constant ratio 
and if the system efficiencies do not change then Sj will have a constant 
value which can be determined. 
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The atmospheric factor 



7.6 


conttibatloo. from oxygen and water vapor tranaialaalvltlee only. 

The operating frequencloa nuat he ehoaen ao that thia la dependent on preaaure 
bat independent of the temperature profile end the water vapor profile. De- 
talla of the aelectlon procedure to aetiefy theaa condltlona are given In 

Section 7e3e 


The background factor 

S * )vi 

characterleea the effecta of absorption by other atnoopherlc constituents and 

sea surface reflectivity, although these by nature vary widely It will be 

cnown in Section 7.2 foUowlng that with reasonable aasuniptlons the product 

s can be made constant when the frequencies are simply related. 

B 

When the conditions established here are met by the appropriate choice of 
operating frequencies and aystem design Chen the Index. S, will be the desire 
measure of pressure. Measurements of other atmoapherlc parameters are not 

needed. 


7.2 Sensitivity to Background 

The assertion in Section 3.3 that the frequency dependence of the re- 
fleeted elgnal strength, due to ocean reflectivity end 
other than by oxygen and water vapor, can be expressed as 1 
be Justified later. However It will be shown here that undesirable back- 
ground sensitivity of this for. can be eliminated from the ultimate deter- 
mination of pressure by combining the measurmi reflected Intenalties at 
appropriately chosen frequencies. This removes the necessity for simultaneous 
subsidiary meaaurementa to determine the unknown and variable coefficients 

a. b and c. 


a 
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Using 


; 0 X 

( 'X " )vj 


^Q(a+bv+cv2) 


7.8 


tt is readily shown by substitution that Sg is independent of the coefficients 
a, b and c which describe the background variation of the following conditions 
on r^^ and are satisfied 


V 

T 


Vr V 
2 i i 

i 


0 


0 


7.9 


Vr V 2 

2 i i 

i 


0 . 


Adopting the ratio measuring technique, where the indices are related by 
equation 7.5, Immediately satisfies the first condition in equation 7.9. 

To satisfy the other conditions a minimum of three pairs of frequencies are 
required in which case the conditions simplify to 


(V2“'’l) + i^(V4~V3) + P(Vf,-V5) = 0 

and 7 . 10 

(v^-v^-) + rtCv^-v-^) + = 0 


where -r 3 = ri, = « and -r^ = rg ® while making -ri = = 1 is not an 

effective limitation on either system design or performance. 


The relationships 7.10 between ix, 3 and the frequencies are purely 
mathematical with an infinite number of solutions. It should be noted that 
the conditions can still be fulfilled if Vj. = V3 and/or Vi, ® V5, so that 
the instrument can have 4, 5 or 6 operating frequencies. 


7.3 Temperature and Water Vapor Sensitivity 

The contribution to tlio strengtli of the return signal, due to oxygen 
and water x-apor absorption has been written In terms of the total integrated 
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absQvption coefficient as Then the contribution of these 

constituents to the atmospheric factor, S^, of the pressure index, S, 
expressed as a product in equation 7.6 can now be put in terms of the sum 
of the coefficients es 

s, . l0-»-2 S|'l''l<''t)| ’-ll 

A i 

The sensitivities of S then to changes in surface pressure, temperature and 
water vapor content (characterlziug profile changes) are thus simply related 
to the corresponding sensitivities of the summation term. 


For S to be independent of the temperature profile we require that 
d 


L i J 


OTg/Tg) 

which becomes, using equation 6.4 

^ri <yi+l> ^''l^ 


- 0 


7.12 


7.13 


Similarly for S to be independent of water vapor content we obtain the 
condition that 

2'i K<v]«, ■ 0 

i 

Equations 7.13 and 7.14 represent further constraints on the choice of 
operating frequencies which, in contrast to conditions 7.10, are dependent on 
atmospheric absorption characteristics. The derivative data of Fig. 7 are 
necessary to enable a set of frequencies related by equations 7.10 simultane- 
ously to satisfy these conditions. 


7.4 Pressure Sensitivity and Signal to Noise Ratio 

Since pressure is the parameter to be determined it is important that 
the sensitivity of the index S to pressure changes should be maximized. The 
fractional change in S for a fractional change 6 in the surface pressure can 
be obtained from equations 7.11 and 6.3 thus 
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^ = -0.2 In(lO) 

It is apparent from an examination of Fig. 7 that this is greatest when the 
main pressure sensing channels are well into the oxygen absorption band, 
■iiowever the very heavy atmospheric absorption reduces the strength of the 
return signal and a consideration of the signal to noise characteristics 
becomes essential. Obviously a compromise between increasing pressure 
dependence and decreasing signal to noise must be achieved. This affects 
the choice of operating frequencies in a way which will depend on the opti' 
mization criterion adopted. 


The potential instrument performance can be analysed with sufficient 
accuracy for the present purpose by considering the noise generated within 
the receiver and omitting the much lower noise contributed by atmospheric 
emission and reflection. It is assumed that the bandwidth of the signal, 

Af , can be adjusted to optimize the system performance and that the detector 
and signal bandwidths are equal. When the received signal is measured 
with a square law detector followed by an integration of time constant t, 
the noise-to-signal ratio at a frequency is 


Ni 




1 

(tAf)'^ 


!i. 


AkTAfF I 

1 


7.16 


and the overall noise-to-signal ratio of a multifrequency instrument is ob- 
tained by adding the squares of the individual ratios thus: 




1 AkTAfF) 

I VVI 


7.17 


If the bandwidth is chosen to minimize this then the minimum occurs 


when 


(AkTAfF)2 



2,| i R i ) 
i 


which by substituting for Af in equation 7.17 gives 


L'T 


T < i 

Assuming that in the absence of atmospheric attenuation the received 
power is the same. in all channels then 

P (v ) » P 10”°*^^T^'’i^ 

i' o 

and equation 7.18 may be written in terms of 

The everell accuracy o. the laetcu«nt can be conveniently d.acrlbed 
by an Index of petfotnance. I. obtained by dividing eguatlon 7.15 eguered 

by equation 7*18 to give 

, (.21n(10))2 52 t Pq . M 


8kTF 


where 


M 


iv* '2 
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7.21 


The Index I Is the ratio of (the change In pressure Index for a fractional 
change In pressure, S) to (noise). The mlnia«.m value of I - 1 determines the 
limiting sensitivity of the Instrument. M is a merit function which depends 
only on the choice of operating frequencies and Che properties of the 
atmosphere. The numerator is a measure of sensitivity to pressure end the 
denominator Is a measure of atmospheric absorption. The system performance 
Is optlmlaed by choosing frequencies which maximise M. 

Another useful Indicator of the sultablUty of a chosen frequency set 
can be obtained from the change In Index S lor a fractional change In the 
mirface pressure. This is the fractional accuracy which the signal pro- 
cessing of the Instrument must be designed to provide. The total Irsctlona 
accuracy will be divided among the channels and When squally distributed 
between n channels this may be written as 


F - 0.2 In(lO) tA^ I 2'i’‘i‘'t^''i’ I 

^ i 


7.22 


-39- 


The value of is determined by the eljoice of frequencies and has important 
consequoncies for system design whicli will be considered in Section 11 when 
the errors associated with sea surface statistics are analysed. Meanwhile, 
it is important to recognize that a relatively large value of will pose 
less technical problems than a small value. 

7.5 Operating Frequencies 

The conditions established in the preceding sections and the deriva- 
tive data of Fig. 7 enable a systematic search to be made for suitable seta 
of operating frequencies. Such a search has been made both numerically and 
by hand with different restrictions on the number and range of operating fre- 
quencies. The results are presented in Table 2 to enable easy comparison of 

the alternatives available. 

Numerical optimization for the selection of up to six frequencies plus 
the two indices, o and $, has been performed with the frequency range re- 
stricted to between 40 and 60 GHz. Over this relatively narrow range the 
smooth frequency dependence of water vapor absorption is adequately repre- 
sented by the background expression and it is not necessary to impose sepa- 
rately the water vapor independence constraint, equation 7.14. The results 
of these calculations. Table 2, line 1, show that the best performance is ob- 
tained using five frequencies. The maximum of 51.2 GHz is effectively limited 
by the requirement for temperature independence. This keeps M low and leads to 
the need for very precise measurement as indicated by the fractional accuracy. 
Fa = 0.05%. When the requirement for temperature Independence is relaxed only 
four frequencies are needed. Table 2, line 2. The maximum frequency at 53.2 is 
further into the band and consequently the merit function is larger. In this 
case, the required fractional accuracy of 0.31% can more easily be met but a 
subsidiary measurement of the atmospheric temperature to an accuracy of 1.2® is 
required to allow the pressure to be accurately determined. 

An indication of the typical behavior of the merit function as the fre- 
quencies are changed is shown in the example of Fig. 10. This contour plot of 
log M for continuous changes in the lowest frequencies and the spread of fre- 
quencies shows that M falls rapidly as the highest frequency increases above 
about 53 GHz. The maximum value of M lies on something of a plateau so that a 


Fractional 



Table 2: Characteristics of Various Sets of Operating Frequencies 



























reasonably wide* ranRo of .sotn of operating frequencies arc available without 
excessive reduction in M. Superimposed on the contour plot la a line of 
temperature Independence and it ia seen that within the range of frequencies 
on the plot maximum M and temperature insensitivity cannot simultaneously be 
obtained and a compromise choice must be made. This plot was obtained using 
smooth derivative data so that detailed structure associated with the indivi- 
dual lines of the oxygen band is not Included. The smoothing was done to 
ensure that the optimization routine did not select a local maximum of M but 
naturally any final choice of operating frequencies would allow for this 
structure. The plot shown here is representative of several which were com- 
puted for different numbers of and arrangements of frequencies and is a use- 
ful Indicator of where to look for appropriate sets of frequencies. It Is ap- 
parent that for maximizing M the spread for frequencies should be as wide as 
possible and the main pressure sensing channel should not be further into the 
band than about 53 GHz. There would also seem to be advantage to be gained 
by extending the lower frequency limit below 40 GHz. 

With the lower boundary set to 22 GHz it is necessary to introduce 
specifically the water vapor Independence constraint equation 7.14, because 
of the effect of the spectroscopic line at 22 GHz. Operating frequencies were 
sought by hand and near optimum sets could be obtained in which the fine 
structure of the derivative data was also used. In this case the requirements 
for temperature and water vapor insensitivities were found to be competing 
constraints. Consequently compromise solutions were sought in which temperature 
and water vapor sensitivities were reduced to reasonable levels. The results 
of this analysis are shown on line 3 of Table 2 and the required knowledge of 
temperature (better than ± 20°C) and total water vapor content (better than 
0.5 gm cm“^) are well within the capabilities of current Instrumentation. How- 
ever, so long as the operating frequencies are restricted to the lower side of 
the oxygen absorption band little improvement in Pressure Sounder performance 
can be expected. 

When the possibility of using one pair of frequencies in the 60 to 100 
GHz range was explored systems with significantly better potential performance 
were obtained. In the six frequency schemes shown in lines 4 and 5 of Table 2 
sensitivity to the temperature profile has been reduced to an insignificant level. 





The water vapor content change of 25 gm cm*"^ which gives a change In the pres- 
sure Index equivalent to a 3mb presnure change la approximately three tlmen 
larger than the 7 gm cm"^ maximum value of this parameter and If necessary 
a correction could be applied using climatological models. The sensitivity 
to pressure as indicated by the merit function and fractional accuracy is 
better than when only the lower wing is used. This is the result of using 
two channels with high sensitivity to pressure. The consequently less ex- 
acting demands on signal processing accuracy allow more flexibility in system 
design as will be discussed in Section 11. 

7 . 6 Summary 

The procedure which has been used for the selection of sets of operating 
frequencies is summarized in the block diagram of Figure 11. The examples 
of Section 7.5 show that sets can be chosen which reduce or effectively elimi- 
nate dependence of the signal on background effects, temperature profile 
and water vapor content. The following section is devoted to the validation 
of the selection procedure which is Indicated at the bottom of Figure 11. 

The selection of an optimum set of operating frequencies is the subject 
of continuing investigation. The numerical optimization routine is being 
extended to include a reformulation of the background effects. No great im- 
provement in the sensitivity to pressure is to be expected since this is 
limited by how close the main pressure sensing channels are to the band 
center. So long as maximum global coverage is required with transmitter 
powers of a few watts, the limit has probably been reached in the final 
examples of Table 2. However gains may be possible in the treatment of tem- 
perature, water vapor and background effects by using slightly different fre- 
quencies or perhaps by using an arrangement of frequencies not yet examined 
although it is not anticipated that more than six frequencies will be used. 

The ultimate decision on which of the possible sets should be. imple- 
mented will depend on other factors which are considered In later sections of 
this report. The advantages apparent in Table 2 which are gained by using 
a wide frequency range must be carefully weighed against the disadvantages 
of a degree of uncertainty in the applicability of the background model over 
the full range and more difficult Instrument design, particularly with regard 
to absolute calibration. 
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Figure 11. Summary of procedure for selecting operating frequencies 
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8 . Numerical Slwulatlona 

8.1 The Use of Radiosonde Data 

The method was detailed In Section 7.3 by which the operating frequencies 
of the pressure sounder could bo chosen such that the measurement of pressure 
was independent of variations In atmospheric temperature and water vapor pro- 
files. The very simple form of the standard atmosphere and the simplest of 
profile variations were used as a basis for the necessary calculations, but 
In practice the atmosphere varies greatly from this stylized profile, parti- 
cularly so for the water vapor content. It is necessary therefore to assess 
the validity of the design procedure by examining whether the desired inde- 
pendences are achieved with a much wider range of profiles. 

Radiosonde measurements provide a data set Ideally suited to this pur- 
pose. The standard radiosonde observations made and recorded on punched 
cards are pressure, temperature, relative humidity and wind velocity for the 
surface and height, temperature, relative humidity and wind velocity at 35 
standard pressure levels from lOOOmb to 2mb. Surface pressure Is recorded to 
the nearest mb, temperature to an accuracy of 0.1°C, relative humditly to 1% 
and height to the nearest 10m. The wind velocity data is not needed for the 
numerical simulations. The standard pressure levels are at intervals of 50 mb 
from 1000 to 200 mb , at 25 mb intervals to 100 mb , at 80 mb , at 10 mb Intervals 
to 30 mb and then for the pressures 25, 20, 15, 10, 7, 5, 4, 3 and 2 mb. 

A test set of radiosonde data from two weather stations was used. The 
data were from a) Cold Bay, Alaska where the surface temperature varied from 
4 C to 18 C and the total water vapor content varied between 0.89 and 4.02 
gm cm ^ and b) Balboa, Panama with corresponding variations of 23°C to 31°C 
and *4.13 to 5.61 gm cm In Figure 12 some of these profiles are plotted to 
Indicate the typical variations. The U.S. Standard Atmosphere is included 
for comparison. 

8.2 Calculations 


Numerical simulation of the experiment was carried out by the following 
method. At each of the frequencies of a set chosen by the procedures of 
Section 7 the transmissivity was calculated for a vertical path through the 








atmosphere defined by a radiosonde observation data set. The absorption due 
to oxygen and water vapor for each layer between the standard pressure levels 
was determined from the spectroscopic models of Sections 4 and 5. The thick- 
nesses of the layers were obtained directly from the height data and cor- 
responding average values of temperature, pressure and water vapor mass 
mixing ratio were found from the radiosonde measurements. The average values 
used were simply the linear means of the parameter values at the standard 
pressure levels bounding each layer. This linear averaging procedure is 
sufficiently accurate since a large number of layers are involved and it was 
the method used when the derivatives for Section 6 were computed. 

An estimate was also included of the absorption by the atmosphere re- 
maining above the highest pressure level for which data was recorded. The 
chosen frequencies are away from the line centers so that this contribution 
is small. It was approximated by the absorption of a layer of thickness 
equal to twice the scale height at 2mb with a pressure of half the minimum re- 
corded level. 

The integrated transmissivities thus computed were combined according to 
equation 7.11: 



The frequency set is one for which the background factor, S„, is 1 and if a 

D 

design is assumed such that the Instrumental factor is unity then S. is also 

A 

the index of pressure, S. 

This procedure was repeated for the large number of radiosonde obser- 
vations of the test set so that a scatter diagram could be plotted of the 
measured surface pressure against the simulated pressure index, S. 

However, the radios*. adc measurement of surface pressure is not the best 
value for an atmosphere characterized by the parameter averages which were 
used to compute the transmissivities. Furthermore, it is only recorded to 
an accuracy of ±0.5 mb and this was found to be comparable to the standard 
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deviation in. the scatter of points. Consequently, the surface 
computed by Integrating through the vertical column sampled by 
The total pressure Is the sura of contributions Ap from each of 
where 


pressure was 
the radiosonde, 
the layers 


Ap 


R T* Pave 
ave 


8.2 


Here g is the acceleration of gravity, R Is the gas constant for dry air. Ah 
the thickness of the layer, p is the average pressure and T* , the average 

Sv C 3 V G 

virtual temperature is given in terms of the average water vapor mass mixing 
ratio 8^^^ and the average temperature, T^^^ as 


It 

ave 


1 + 1.609 8 


ave 


1 + 8 . 


ave 


ave 


8.3 


The factor of 1.609 derives from the ratio of the mean molecular weight of dry 
air and the molecular weight of water vapor. 


In Figure 13 the measured surface pressure is plotted against the com- 
puted value. The systematic difference of about 1.2mb is attributable to the 
inadequacies of the parameter averaging procedure. The distribution of 
points about the best fit line has a standard deviation of about 0.4mb and 
most of this is probably rounding error associated with the integral repre- 
sentation of the measured surface pressure. 

8.3 Results 

The six frequency system defined in line A of Table 2 was evaluated in 
this way and the resulting scatter diagram of measured surface pressure against 
simulated pressure index, S is plotted in Figure lA. This diagram shows that 
the relationship between the index and surface pressure is satisfactorily re- 
presented by a straight line over the range of pressure from 950 to 1020 mb. 

All of the points lie within il.Smb about this straight line so that the pres- 
sures as deduced from the simulated signals are within ±lmb of the actual 
surface pressures. From the Gaussian distribution curve fitted to the histo- 
gram of errors, inset in the figure, the standard deviation of the deduced 
pressure can be estimated to be about 0.6b mb. In Figure 15 the error in 
pressure (deduced minus actual to the nearest mb) is plotted as a function of 
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Surface Pressure measured by radiosonde, mb 


Figure 1!^* Computed surface 


pressure vs. reported surface 


pressure. 







Figure 14. Scatter Diagram. Numerical Simulation of Index From a 6 Frequency 3 Ratio 








both ourfaco tamperatuta (rig. 15a) and total watar vapor rontant (Fig. 15b). 

It can ba aaan that no ootralatlon axiata and that the daalgnad Indopandanaa 

o the maaautaniant from tamparatura and watar vapor appears to have bean 
achieved . 


Whan the aomputad rather than the measured surface pressure Is plotted 
on a scatter diagram against the simulated Index, Figure 16, the scatter of 
points is noticeably reduced. This reduction In scatter reveals an offset 
of the Panama data points which Is Indicative of a small residual tempera- 
ture and/or water vapor sensitivity. If due solely to temperature effects 
t en the sensitivity is such that a 30° change In surface temperature will 
change the pressure Index by an amount equivalent to a 3mb pressure change 
while If solely a water vapor effect then a 7 gm crn'I change Is equivalent 

to a 3 mb change. For these levels of dependence a correction with cllma- 
tological models is satifactory. 


However, with a slight modification to the operating frequencies the 
required Insensitivities can be restored. Using the six frequencies of 
line 5, Table 2. the pressure Index was computed for the test set of radio- 
sonde data. The resulting scatter diagram. Figure 17, exhlhlts no tempera- 
ture or water vapor sensitivity and the standard deviation In the deduced 
pressure Is about 0.4mb. Hhile some of this scatter Is probably a computa- 
tional effect due to summation over a finite number of layers there Is evi- 
dence that variability In the atmospheric structure Is also a contributory 
factor. For example, the point most seriously In error, marked A on Figure 

17, corresponds to the profile A In Fig. 12 which has a temperature In- 
version at 7 km. 


These numerical simulations employing a wide range of atmospheric 
structures have demonstrated that frequency sets may be chosen so that a 
pressure index can be determined which Is sufficiently Insensitive to the 
temperature and water vapor profiles that subsidiary measurements of other 
meteorological parameters are unnecessary. It has also been demonstrated 
that a measurement error of not more than 0.4mb, attributable to variability 
of the atmospheric structure which cannot be modelled, can be achieved with 
an appropriate choice of operating frequencies. Refinement of the selection 
procedure may enable this error to be decreased further while sophistication 
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of the computationnl analysis may reduce its upper limit. This atmospheric 
variability error is one contribution to the total error in surface pres- 
sure as measured by the microwave sounder. Other sources of error will be 
considered in Inter sections of this report. 


9. 


Tlie Effecta of Other Atmospheric Constituents 


The presence in the atmosphere of absorbing gases, other than those 
already considered, has a slight effect on the choice of operating fre- 
quencies. In contrast, atmospheric aerosols have considerable Influence 
on the optimum frequencies as well as on other aspects of the system design. 

Several minor constituents have spectral lines in the 25 to 75 GHz 
band of Interest (Waters, 1976) but only ozone has absorption through a 
zenith path of sufficient strength that it may possibly affect the measure- 
ment. In addition some attention must be given to the spectrum of the 
0 I 60 I 8 isotopic form of molecular oxygen so that the line centers may be 

avoided . 

Cloud and rain are the principal aerosols and of these non-precipi- 
tating cloud is of greater concern because of its more frequent occurrence. 

The particular advantage that a microwave sounder has over an optical in- 
strument is its ability to penetrate clouds. Nevertheless, clouds do 
absorb millimeter waves and both the degree of absorption and its variation 
with frequency affect the system design. The assumption made in Section 3.3 
concerning the form of the frequency variation will be justified and the 
absolute value of absorption, which is needed to determine the signal-to- 
noise errors associated with a variety of cloud cover conditions, will be 
obtained from an analysis of scattering by aerosols. 

9 . 1 Ozone 

The mixing ratio of ozone in the atmosphere varies with altitude and 
exhibits a maximum at between 20 and 30 km. Although the total zenith 
path attenuation is small, frequencies near the line centers should be 
avoided because the ozone content and hence the absorption varies by up 
to 50% daily, seasonally and spatially. The strongest ozone lines between 
20 and 75 GHz are listed in Table 3 (Depannemaecker et al, 1977) with their 
peak absorption coefficients. The strongest line is at 67.356 GHz and at 
this frequency the Integrated ozone zenith opacity is typically .04db based 
on an Integrated depth of 20 cm. By comparison oxygen zenith opacity is 
nearly 5db at the same frequency. This is the only line of any real concern 


Frequency, 

Absorption, 


nepers km ‘ 

GHz 

at 220*K 

23.860 

2.3 

28.960 

2.1 

30.052 

4.6 

30.181 

2.6 

36.022 

3.4 

37.832 

7.2 

42.832 

3.9 

43.653 

6.8 

44.871 

2.1 

50.034 

2.6 

51.976 

6.3 

53.688 

14.3 

55.356 

6.0 

58.094 

5.1 

61.347 

4.5 

61.927 

16.4 

63.072 

2.7 

65.236 

16.9 

66.059 

5.9 

67.250 

4.0 

67.356 

48.4 

68.421 

3.1 


Table 3. The strongest ozone lines in 20 to 75 GHz 
range with approximate peak absorptions at 
pressure of one atmosphere. (From Depannemaecker 
et al., 1977 and Kakar, R.K. , private communi- 
cation. ] 
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A Lorentz line shape la generally used and the llnewidth parameter 
has a value of about 2.3 MHz mb'^ (Waters. 1976). Calculations of the 
zenith attenuation are shown in Figure 18 for various ozone profiles. At 
frequencies greater than 100 MHz from the line center the absorption is 
reduced to leas than 1/4 of the peak value. For the operating frequencies 
in Table 2 the closest are 80 MHz and 150 MHz from the line center. Ozone 
absorption for these is respectively about O.OlOdb and O.OOSdb while a 3mb 
pressure change causes a change in attenuation of about 0.03db. A 
correction to the pressure measurement based on the known global distri 
bution of ozone would probably reduce ozone associated errors by a half. 

In the first case this is still significant at O.Smb but in the latter case 
the final frequency set of Table 2, the error is reduced to an acceptable 

level of about 0.2mb. 


9.2 


The isotopic species of oxygeh. end “o. are present In the 

atmosphere in the following fixed 

relative abundancies: 

I60 

99.76 % 

17o 

.037% 

I80 

.204%. 

Thus the relative abundancies of 

the principal molecular forms of oxygen 

are 


I602 

99.52% 

16Q10O 

.41% 

16ol7o 

.07% 


Each of these has a distinct spectrum but in Section 5 and 6 the spectrum 
of only the dominant ^®02 molecule was used as the basis for the frequency 
selection calculations. Since absorption accuracies at the level of 
fractions of 1% are needed to satifactorily measure pressure the proportion 
of I60I80 present in the atmosphere is possibly sufficient to affect the 

calculations. 


% 


Zenith attenuation, db 



Thfi spectrum of differs from th<it for ^®0j> in huvluR twice the 

number of lines. This is because the lack of symmetry allows for even as 
well ns odd values for the rotational quantum number N. As n consequence 
the line intensities for are approximately half as ntrong as those 

of the correspondinR transitions for Accurate values for the reso- 
nance frequencies of ’.iave been published for only a few lines near 

the band center (Stelnbach and Gordy, 1975). Typically the lines for 
are displaced from the corresponding '^0^ lines by a few hundred MUz. Since 
there are also lines for even N resonance peaks will bo scattered throughout 
the ^®02 valleys. However, in the band wings, where the. oporntlng Cre- 
quf... ;<38 of the preasure sounder have been chosen, even the *^02 pouku are 
not very prominent. It can therefore be expected that little if any effect 
of the isotope substitution will be noticeable except perhaps within a few 
MHz of the line centers. Calculations are being undertaken using 

the molecular constants of 0 to determine the positions of other line 

centers so that these may be avoided if absorption effects are found to be 
significant. 

9.3 Absorption by Clouds 

9.3.1 Water Clouds 

The presence of cloud in the vertical path through the atmosphere will 
affect the return signal through absorption and scattering of the Incident 
radiation. Ikich theoretical work has been published using electromagnetic 
scattering theory to compute attenuation coefficients, y, and backseat ter Ing 
cross sections, a (Ryde and Ryde, 1945, Goldstein, 1951, Gunn and East, 1954). 
Rigorous experimental verification of the following theoretical calculations 
for cloud absorption (as well as other aerosol effects) is generally un- 
satisfactory because it is not possible with a reasonable number of measure- 
ments to fully specify the spatial and temporal fluctuations in aerosol 
density and drop size distribution. However, the measure of agreement which 
is obtained, is commonly accepted as evidence that the theory provides an 
adequate basis for the description of electromagnetic wave Interaction with 
aerosols. 

For a theoretical analysis of the absorption by clouds at the fre- 
quencies of interest here, the Rayleigh scattering approximation is valid 
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sl’ice the wavelengths are much larger than the diameter of the water drop- 
lets in cloud (<0.05 mm). Frequency dependences in y and a arise from this 
and also from the dielectric properties of the scattering water droplets. 

To complete the quantitative description of absorption by clouds, assumptions 
must be made about the distribution of drop sizes. Laws and Parsons (1943) 
made extensive measurements and their results are most commonly used. 


With an analysis using these fundamentals Gunn and East showed the 
attenuation coefficient, y, due to water cloud to be given, in db km"^ by 
the expression 


Y = 0.434 ^ Im(-K) 

’ A p 

where the wavelength, X, is in cm, the density of water, p, is in gm cm ^ 
and the density of water in the cloud M is in gm m There is a frequency 
dependence within Im(-K) which is given in terms of the real and imaginary 
parts, e’ and e", of the complex dielec^ric constant, [e], as 


Im(-K) 


3e*' 

(e' + 2)2 + e"2 


9.2 


and details are given in Section 10 of the Debye model from which e* and e 
can be obtained. 

Benoit (1968) fitted straight lines through Gunn and East’s discrete 
data shown in Figure 19 and arrived at the following empirical expression 
for the absorption coefficient; 

Y • e^^ db km”^ 9*3 

The frequency, v. is in GHz, the frequency index, bi, for water cloud has a 
value of 1.95 and the temperature coefficient, a\ “ -6.866 (1 + 0.0045T) 
when T is in °C. More recent calculations by Crane (1971) extend the data 
available to 100 GHz, Figure 20. 

In Section 3.3 it was assumed that integrated cloud absorption could 
be modelled by a second order polynomial variation with frequency, viz 

r. • a + bv + cv2, 

A 
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Figure 19. Attenuation coefficient for clouds 
(From Gunn and East, 1954) 
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Figure 20. Attenuation coefficient as a function 
of frequency for liquid scatterers 
(From Crane, 1971) 






and It can readily be shown now that the above data fits such an expression. 

In Table 4 values for the absorption coefficient, calculated assuming the 
empirical relationship, 9.3 above, are compared with values from a best fit 
polynomial. The values agree to better than 0.2%. Cranes smooth curves in 
Figure 20 can also be fitted, within the 25 to 75 GHz range by other second 
order polynomials and the computed lines are Indistinguishable from the 
original lines. Of course, the actual values of a, b and c are not important 
since selecting frequencies which satisfy the conditions 7.9 makes the pressure 
index, S, independent of these coefficients. 

The magnitudes of the absorption coefficients for various cloud types 
are important since these will determine the strengths of the received 
signals. Thus heavier clouds will produce a reduction in signal-to-nolse 
ratio which ultimately limits the accuracy of the pressure measurement. 

Several models of water clouds are commonly used, in which the cloud is 
characterized by the density, M, and its thickness. The product of these 
two factors gives the Integrated liquid water content, Wj^, usually expressed 
in gm cm~2 and in accordance with equation 9.3 attentuaticn is directly 
proportional to Wj^ at any fixed frequency. Cloud temperature is also 
significant, being generally in the range -10°C to +20°C and varying with 
cloud altitude, latitude and season. Table 5 shows attenuation for a re- 
presentative set of cloud types at 0°C for the frequencies 25, 50 and 75 GHz. 
The values will be slightly less for higher temperatures. 

For more than half the time, atmospheric liquid water is essentially 
zero (Staelln et al, 1976) while the overwhelmingly dominant cloud types 
have Wj^ < 0.05 gm cm Occasionally liquid water exceeds 0.6 gm cm“^, the 
highest example in the table but these are generally small scale phenomena. 

The presently envisaged pressure sounder designs should allow meaningful 
pressure measurements to be made with attenuations of up to 5db in the 
highes*- frequency channel. Since this corresponds to Wj^ 'v. .1 gm cm”^ almost 
all cloud conditions will be accommodated. Further discussion of potential 
coverage is reserved for Section 12 following a complete description of the 
instrument and an analysis of signal-to-nolse errors in Section 11. 


Frequency 


Cloud Absorption, db 


From 

Y*(constant)xv 


From Y" 

-4.412 X 10“2 
+4.175 X 10“^v 
+8.412 X 10"S2 


0.5869 

1.1312 

1.8466 

2.7309 

3.7825 

5.0000 


0.5860 

1.1325 

1.8472 

2.7302 

3.7814 

5.0008 


Table 4. Approximation of cloud absorption by a second order 
polynomial • 


Density 
gm m”^ 

Depth 

km 

Integrated 
Liquid Water 
gm cm“^ 

Absorption, db 



25 GHz 

50 GHz 

75 GHz 

Stratus C 
0.1 

!louds 

1.0 

0.01 

0.05 

0.21 

0.48 

0.25 

1.0 

0.025 

0.13 

0.54 

1.21 

Fair Weather Cumulu 
0.5 1.0 

B 

0.05 

0.27 

1.07 

2.41 

1.0 

1.0 

0.1 

0.54 

2.14 

4.82 

Tall Cumulus 
1.0 2.0 

0.2 

1.07 

4.29 

9.64 

2.0 

3.0 

0.6 

3.22 

12.9 

28.9 

Cirrus 






0.05 

0.5 

0.0025 

.00006 

.00012 

.00018 

0.10 

1.0 

0.01 

.00023 

.00046 

.00070 


Table 5« Integrated absorption for some typical clouds 
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9.3.2 Ice Clouds 

High altitude cirrus clouds frequently occur and will often extend 
for 1000 km or more. These clouds consist of small ice crystals, are of 
low density, M < 0.1 gm m"^, and are generally thin layers, 1 km. The 
dielectric constant for ice is very much different from that for water 
and this results in a much lower absorption coefficient. Gunn and East's 
data for ice cloud attenuation are also shown in Figure 19 and Benoit 
fitted these with an empirical expression of the same form as equation 9.3. 

The indices for ice clouds have the values bi * 1.006 and 

ai • 8.261(1 - 1.767 x10”2t- 4.374x10 "*' t2) 9.5 

This is very close to a linear vari*»tion with frequency and is consequently 
readily accounted for by a second order polynomial. However, the values 
in Table 5 for attenuation by typical ice clouds, at 3 or 4 orders of 
magnitude less than for water clouds, Indicate that the presence of cirrus 
will have an insignificant effect on pressure measurement. 

9.4 Absorption by Precipitation 

The distribution of drop sizes in rain typically exhibits a maximum 
number of drops with radii near 0.5 mm (Marshall and Palmer, 1948) and, 
depending on the rain rate, some drops may have radii up to 4mm. The 
wavelengths most suitable for pressure sounding are 'v< 4 mm and longer. 

Because of the similarity of these dimensions it is necessary to use Mle 
scattering theory to determine the absorption coefficients for rain. The 
results of such calculations by Crane are included in Figure 20 and the 
smooth cusrves can be adequately fitted with a second order polynomial in 
frequency. Absorption by rain can therefore be accounted for provided the 
total attentuation is not excessive. Typical rain cloud states must next 
be examined to establish the conditions under which this requirement may be 
satisfied. 

Figure 21, taken from Valley (1965) shows how the liquid water content 
of clouds varies with altitude for surface rain-rates between 2 and 20mm hr"'. 
The total atmospheric liquid water in these cases can be estimated to vary 
from 0.04 to 0.15 gm cm"^, partly in the form of small cloud aerosols and the 
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remainder as larger rain drops. Without more detailed information on the 
vertical extent of the rain it is impossible to evaluate with any accuracy 
the total attenuation but assuming, for the 2mm rain rate, a 3 km rain 
depth the total rain and cloud attenuation can be estimated using Figure 20 
to be less than 5db at 75 GHz. Other precipitation conditions occur with 
rather thinner cloud of perhaps higher density which could result In similar 
total absorption. It is possible then that the pressure sounder may give 
meaningful measurements with rain rates of a few millimeters per hour. Its 
use in heavier rain appears to be limited by attenuation rather than in- 
accuracies in the frequency variation model. 

More detailed and accurate analysis of absorption by rain is not neces- 
sary at this time since an instrument whose use was limited to non-preci- 
pitating clouds would nevertheless be worthwhile. Similarly, although models 
exist for absorption by other forms of precipitation which would allow a 
superficial assessment, a detailed consideration of the possibilities for 
measurements in hail, sleet and snow is of no great value because of their 
comparatively infrequent occurrence. 

9.5 Backscatter from Aerosols 

If it is sufficiently strong, backscatter from aerosols could affect 
the measurement of surface pressure by this proposed microwave transmisso- 
meter technique since the return signals would not have been attenuated by 
the Important lower levels of the atmosphere. However, the following analysis 
shows that cloud reflectivity is not significant for those conditions when 
attenuation is small enough to allow the pressure measurement to be made. 

The reflectivity of rain has been used extensively in weather radars to 
monitor the position and extent of precipitation. At millimeter wavelengths 
the backscatter cross-section is large enough to affect the pressure sounder 
and will probably be a more severe limitation on measurement capability than 
attenuation by rain. 

The electromagnetic scattering problem was solved using the Rayleigh ap- 
proximation by Gunn and East (195A) to give the magnitude of the radar 
signal reflected by cloud as 

Pr - Pt a L |k|2 . 

8 h2x‘* 




9.6 
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Most of the symbols have been defined in Section 3. L Is the depth of the 
cloud, D Is the particle diameter and the summation Is over the distribution 
of drops In unit volume. There Is an Implicit frequency dependence In |kP 
which Is given in terms of the complex refractive index [e] by 


. tel - 

te] + 2 • 


9.7 


It Is convenient for our purposes to define a normalized cloud backscatter 

cross-section, o , which may be directly compared with the sea surface cross- 
* o 

section, as used in equation 3.4. At the same time we follow convention 
by Introducing the radar reflectivity factor 


Z 



9.8 


Then 


" ttS IkP L Z X 10“^ 
X** 


9.9 


and Is a dimensionless ratio when a confusing array of units Is used: L In 

km, Z by convention In mm^ m~^ and X In cm. Using the Laws and Parsons 
drop-size distribution for cloud Z can be computed and this results In the 
following empirical relationship with M, the cloud liquid water density In 
gm cm“^ (Valley, 1965) 


Z - 0.048 m2 . 9.10 

Now we are In a position to evaluate a^. Sufficient Indication of the 
possible effects of cloud backscatter can be obtained with an estimate of 
0^ for cloud conditions some five times worse than can be penetrated by the 
pressure sounder. At the highest operating frequency X ■ 0.4 cm and IkP^uO.S, 

depending on temperature. For a deep, heavy, essentially opaqur cloud with 

-3 -<> 

M ■ 1 gm m and L ■ 5 km, equivalent to a total water content ■ 0.5 gm cm 

then equation 9.9 gives 


o - 1.43xl0‘** or -38db. 9.11 

c 

This is Insignificant when compared to the sea surface cross-section, 

0 ° 'v* 8db so that backscatter from water clouds will not affect the measure- 


ments . 
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Little data is available on the reflectivity of cirrus (ice) clouds 
but Valley gives an order of magnitude estimate of Z, viz 

Z 570 M 9.12 

Particularly bad cirrus clouds might have M ■ 0.1 gm m~^ and L “ 1 km 
so that at A « 0.4 cm for which |Kp 0.2 equation 9.9 gives 

- -25db. 9.13 

This is more than three orders of magnitude less than o° and will therefore 
not influence measurements at nadir. Off nadir measurements may occasionally 
be affected hy cirrus cloud reflections if equation 9.12 is accurate. 

For rain, the backscatter cross-section a is also given by equation 
9.9 except that Z must be replaced by Z^, an equivalent radar reflectivity 
factor. This is necessary because the Rayleigh approximation becomes in- 
valid and Z must be computed using the exact Mie equations with the Marshall- 
Palmer rain drop-size distribution. Valley has surveyed the data relating 
Zg to the rainfall rate R (mm hr~^) and at millimeter wavelengths a good 
approximation is 

= 300 R (mm®m“3) 9.14 

Thus, when R ® 5mmhr“^ and L = 1 km the normalized backscatter cross section 
at a wavelength of 0.>'. cm is 


Oj. = 3db. 9.15 

This is consistent with the results from Godard (1970). Being of the same 
order of magnitude at o° it can only be overcome with a pulsed radar and 
precision range-gating. The marginal improvement In coverage which would 
be provided by this technique does not by itself warrent the increased com- 
plexity and expense. 

Sleet, snow and hall also reflect at millimeter wavelengths. For 
these the equivalent radar reflectivity has a similar value to Z^ for rain 
and measuring capability will be similarly limited. 
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9.6 Broken Cloud and Other Non-Uniformities 

All of the design considerations to this point have implicitly assumed 
that the properties of the field of view, smeared for the integration time 
for n single pressure measurement , are constant and hence that the signal 
levels are constant, except for statistical fluctuations due to the sea 
surface. Any non-uniformity along the orbit track such as pressure changes, 
temperature or water vapor profile differences, the presence of broken 
cloud and sea state changes will cause the strength of the return signal 
to vary. Of these the most Important is broken cloud since it can produce 
sudden large changes in the echo signal from nearly adjacent points while 
the other factors will cause smooth changes. When the cloud sizes are 
comparable with the instantaneous field of view then the received signal 
will fluctuate along the track. However, if the cloud sizes are much 
smaller but extend over a large area then the changes will be smooth. In 
either case the effect is to modify the frequency dependence of the back- 
ground as modelled by equation 7.8. An analysis of this problem indicates 
how all non-uniformities may be accounted for since they can all be in- 
corporated into the background effect. 

It has been shown that the appropriate choice of frequencies will 
make the background factor, S„, unity when the cloud absorption coeffl- 
dent is represented by the second order polynomial, equation 7.8. When 
uniformly thick cloud covers a fraction p of the field of view smeared 
along the orbit track then the cloud transmissivity, t^, can be written 

Ta^ (Vi) - (1-p) + p 10 ^ 9.16 

The background factor is obtained from the combination 

h • n |V <>i)| . 9.17 

This will be unity only when p > 0, 1 and can be expected to have the 
largest error when p ■ O.S. 

It was Indicated in Section 9.3 that the worst cloud through which 
measurements would be possible will attenuate the signal by 5db at 75 GHz 
over a one-way path. Using the values of the coefficients a, b and c 
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glven In Table 4 for this cloud Ptnte and the frequenclea of Table 2, line 5 
the calculated value of when p * 0.5 is 

Sg = .9763. 

This represents an error in the measured index S of 2.4% which is equiva- 
lent to a 3mb pressure error. This is a significant error but it 
the worst case and generally errors will be much less or zero. For 
example, when the cloud attenuation is Idb at 75 GHz and p is again 0.5 

the error is < 0.5 mb. 

The effects of broken cloud can be mitigated in several ways. The 
most easily implemented method is to divide the integration time into several 
shorter periods, calculate the index S for each subdivision and finally to 
average S over the total integration time. Processing the measurements in 
this way means that errors only occur when the subdivisions are partially 
filled with cloud. Alternatively, or indeed additionally, fluctuations 
in the return signal strength due to broken cloud might be utilized to 
apply a correction to S. 

Other possibilities involve modifications to the procedure for 
selecting the operating frequencies. A different function of frequency 
to describe the background effects may be advantageous. Water cloud is 
the most Important background effect and it has a nearly quadratic fre- 
quency dependence (equation 9.3), Thus a three parameter model for back- 

ground with the form 

r, “ a + bv^ + cv^ 9.18 

A 

would adequately account for uniform cloud and at the same time signifi- 
cantly reduce the effect of broken cloud. The alternative to this is to 
extend the numerical optimization routine, as indicated in Section 7.6, so 
that it includes a cloud effect term. This can be used with a set of 
cloud states to select operating frequencies which minimize sensitivity to 
uniform and partial cloud cover conditions. 


All of these methods for reducing the broken cloud errors are the 
subject of continuing studies as are the effects of smooth variations in 
the other factors. A quantitative assessment of the relative frequency 
of various atmospheric states is necessary to enable a decision to be 
made on which- of the alternative techniques is most appropriate. 
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10 . S oa Surface. Rnflcctl vlty 

10.1 Normalized BacHscattor Croas-Soctlon 

Rof lection of oloctromaBnetic waven by the Hoa aurfaco In characterized 
by the normalized backscatter crosa-aection, o(0), which ia defined as the 
ratio of (the radar backacatter croaa-aectlon) to (the horizontal surface 
area intercepted by the antenna beam ) (Tomlyaau, 1974). BelnR a ratio. 

0 (G) ia dimensionless and physically represents how good a reflector the 
sea surface is In comparison to a perfectly reflecting sphere. The average 
value of o(0) Is a function of surface roughness, complex dielectric con- 
stant. temperature, frequency, polarization and angle of incidence. Of 
equal importance to the average value of o(0) is the statistical variation 
in the reflected signal due to the nature of the sea surface and this is 
considered in Section 10.8. Many papers have been published in recent 
years on both the theory of sea surface reflectivity and its experimental 
measurement (see Barton,1975, Skolnik, 1970 for bibliographies). 

Measurements of the angular dependence of a(9) confirm the theoreti- 
cal results of Barrick (1968) shown here in Figure 22. With transmitter 
powers limited to a few vatts it will only be possible to make measure- 
ments at or near to the vertical. The nadir value, o°, is greatest for 
low values of the mean sea slope, s. For these calm sea states o(0) falls 
off rapidly as the angle from nadir increases. As the sea becomes in- 
creasingly rough the maximum value Is reduced but o(0) falls less rapidly 
away from the normal. Consequently rough sea states are better for 
making measurements away from the subsatellite track. The rapid fall-off 
of 0(0) at angles greater than 15° for commonly occurring t.ea states 
limits the use of the Microwave Pressure Sounder to nadir and near nadir 

operation. 

10.2 Theory 

In the composite model commonly used the slope of the sea surface 
is characterized by two roughness scales corresponding to gravity and 
capillary waves. Significant waveheights, H, associated with the first 
range typically up to 10m while for the latter the characteristic length 
is of the order of centimeters. The sea surface is considered to be rough 
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Figure 22. Averngo backscatter crona- 
aection per unit area, o“, 
VH Incident angle, a la 
r.m.a. roughness slope, 
Gaussian slope probability. 
(From Barrick, 1968) 


Figure 23. Variation of o® with fre- 
quency for different values 
of 6 (vertical polarization). 
(From Wiltse e_> al, 1957) 


Figure 24. o® as a function of frequency 
(From Grant and Yaplee, 1957) 
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If H > X/2it, where X Ic the wavelength of the Incident radiation, while a 
slightly rough surface has H X/2ir. Both scales are present at the same 
time in the sea surface and the superposition of mlcrostructure on the 
large scale roughness obviously affects the distribution and size of the 
specularly reflecting facets. However, scattering by each scale may be 
analyzed separately and to a first approximation the total backscatter 
cross-section is the sum of the two individual contributions (Bartick and 
Peaks, 1968). At normal incidence the Influence of the larger scale 
dominates since backscatter due to the gravity waves is 3 or 4 orders of 
magnitude greater than that due to the capillary waves. The operating 
frequencies of the microwave pressure sounder are between 25 and 75 GHz 
(wavelengths between 12 and 4 mm) so that even on the calmest days the 
sea surface will be "rough." 

For rough surfaces the recognized technique for the derivation of 
0 ° is the specular point theory (Barrick, 1968). Various formulations 
(e.g. Kodls,1966; Stogryn, 1967) of the problem arrive at the same result 
for the backscatter cross-section which, with the assumption of a Gaussian 
probability distribution of the sea slopes (or heights), Barrlcks gives as: 

a° = |r(0) P 10.1 

s^ 

for normal incidence, where s is the mean sea slope and R(0) is the Fresnel 
reflection coefficient at normal incidence. o°/|R(0)p is independent of 
frequency and Barrick' s values of this lie between 5 and 27 db for various 
sea states while |r( 0)P is between -2.5 and -4.5 db. At normal incidence, 
0 ° does not depend on polarization and the reflected signal suffers no 
depolarization. Cox and Muok (1954) showed that the distribution of wave 
heights is approximately Gaussian and this has most recently been verified 
by Weissman and Johnson (1977). However, Barrick and Snider (1977) show 
that the Gaussian assumption is not necessary while Barrick (1968) shows 

that other statistical distributions lead to similar values and properties 

r o 
of 0 . 


Although the contribution to the total backscatter cross-section 
from the surface mlcrostructure is small a brief consideration of it is 


desirable. In the circumstances where the surface la slightly rough, l.e. 

H < X/2it, 0 ° Is determined by the boundary perturbation method of Rice 
(1951). This leads to a Bragg scattering criterion which make o(6) dependent 
on frequency, angle and polarization. However, at normal incidence o /|r( 0)|^ 
is- again independent of frequency and polarization as the measurements of 
Gulnard et al. (1971) show. 


10.3 Temperature and Frequency Dependence of o° 

The temperature and frequency dependence that there is in a at normal 
incidence arises because of the inclusion of R(0) in equation 10.1. The 
Fresnel reflection coefficient may be expressed in terms of the complex 
dielectric constant, [e ] . 


|R(0)P 




1 + [e]^ 


10.2 


e at microwave frequencies is given by the Debye expression, which, in 
its most general form is: 


e ■ + 

00 


1+ (jWT)^~“ 


0 

Wty 


10.3 


0 ) « 2irv is the radian frequency with v in Hertz 

e is the dielectric constant at infinite frequency (=4.9 ± 20%) 

00 

e_ is the static dielectric constant 

T is the relaxation time in seconds 

a is the ionic conductivity in S m ^ 

a is an empirical constant (best value “ 0.02 ± 0.007) 

e - 8.85 X lO''^^ Fm"! ig the permeability of free space 

0 

Klein and Swift (1977) have reviewed the recent work done to establish 
accurate microwave values for the real and imaginary parts of e demanded by 
advances in precision r.icrowave radlomctry. Uncertainties exist in the 
absolute values at millimeter wavelengths because of the high error asso- 
ciated with e but of principal concern in the design of the pressure sounder 

OO 


ni 




is the variation with frequency expressed explicitly by equation 10.3. 

The quantities c^, T and a are temperature and salinity dependent. 

Klein and Swift (1977) give polynomial expressions from which they may be 
obtained. For the purpose of identifying the frequency dependence of 1 r( 0)12 
from equations 10.2 and 10.3 the following values at 15*^0 and 35%o salinity 
were used: = 74.10. t = 10.50 x lO'^^g and o = 4.290 Sm ^ In addition 

o was approxiLted by zero. The values of IrCO)!^ thus determined for v - 20 
through 80 GHz are given in Table 6. The frequency variation can be adequately 
represented by the second order expression 

+ bv + cv^) 

where a = -.1696, b = -.00294S and c » 4.163 x lO"®, as illustrated by the 
values based on this approximation given for comparison in the third column. 

Temperature dependence is such that at 50 GHz 1 r( 0)| .45/ and .523 

at 5°C and 25°C respectively. Such variation is not significant in the con- 
text of the proposed Instrument where the values of the coefficients a, b 
and c are eliminated by the combination of transmissivity ratios. 

10.4 Measurement of 

Of the many measurements that have been made of sea surface backscatter 
cross-section few have been in the frequency range of interest here. The 
first such measurements reported are those of Wiltse et (1957) taken 

with a 17° beamwidth radar, mounted on the bow of a ship, at frequencies of 
9.6, 24, 35 and 48.7 GHz for a range of angles and sea states. The angular 
dependence has the form given by theory and the measurements at and near 
normal are summarizad in Figure 23. Wiltse ^ conclude that o° does not 
vary with frequency - at least to within the limits set by their 2.5db 
accuracy. The normal incidence values of o'’ at 48.7 GHz is 7 db for a wind 
speed of 4-10 knots. 

Making measurements at 9.4, 24 and 35 GHz from a bridge with beamwidths 
of about 3°, Grant and Yaplec (1957) obtained a normal Incidence significantly 
higher values of o° and (Figure 24) a variation with frequency which has not 
since been corroborated. 


Frequency 

GHz 


1r( 0)P from 
Equation 10.2 


0.5928 

0.5574 

0.5241 

0.4937 

0.4661 

0.4411 

0.4184 


1r(o)|2 

Approx. 


0.5932 

0.5569 

0.5239 

0.4938 

0.4663 

0.4412 

0.4182 


-2.27 

-2.54 

-2.81 

-3.07 

-3.32 

-3.56 

-3.78 


Table 6. Variation of the Fresnel Reflection Coefficient with 
Frequency 


Surface Wind Speed 
(knots) 


0° (db) 


12.5 

25 

29.5 

36 

45 

10 

13 

3 

8 1 

9 


Table 7. Normalized Backscatter Cross-Section at 13.3 GHz 
5° Angle of Incidence (From Krishen, 1971). 




















Uter work ha. conoaotrarod o„ c„„f, rating rha thoorntlrnl angnlnr Jc- 
paodaora of tha baakaaarrar with tha Intantlon of relating thl. to nan „tnta 
para^atara and wind apaad. Gulnard at nl. (1,71) wnda nultlfrannanty 
(0.4 to 9.0 Olln) maaauraaanta bat abaoluta valuaa of o» warn not datamlnad 
^an„ra»anta by Krlahan (1971, with a fan baa. at 13.3 0„r, are given ’ 

.bla 7. Thaaa ara for an angle of 5° fro. nadir but should not ba varv 
arant fro. normal Incidence values. At tha high wind speeds for whirl. 

but ““'ll' l>«ve significant fan. covarnna 

but no indication of this la reported. Barrick (1974) reports measuramant, ’ 

y enest at 9.0 GHr to an accuracy of 1.9 db (Figure 25). It can be 

seen that as the wind speed Increases o” reduces to a mini.... value of about 

the 11 represented by 

et al 0977.° "“r"«"“rs (1.7 db) are those ™de by Jones 

M iiTmal 1 l-5°- Their r.sults 

fl, . . 7 “ =“”-*rlied in Table 8 and In Figure 26. These again 

foam ° "S**" Crnbable 

roam coverage was not reported. 


Measure.ents have also been »de of radar backscatter cross-sections at 
. 3p to 10.6U fro. rough .atalllc surfaces, cheo and Renau (1969) show that 
te specular point theory can be applied at these very short wavelengths and 
at o is independent of frequency for rough surfaces satisfying the con- 
Itlon h/i > 1/4. For a surface with ms height fro. the mean hn,7„ and tha 
mean scale size along the surface ln,50p, n” is 15 db and when hn.lp and l-vlo„ 
then o is 17 db. These correspond to the results of Barrlck-s theory (Figure 22) 
when the surface slopes are s = tan 10° and s - tan 8° respectively. 


10.5 Sea State and Foam 

The roughness of the sea Is driven by eh. wind but the relationship of 
SM state to wind Is rather Imprecise. This is because gravity waves take on 
the order of hours to build up and decay while the capillary waves are ^.ch 
eore Inunedlately responsive to wind variation.. The snu,oth eurve of Figure 27 
from Barrlek (1974) relating sea surface slope to wind speed disguises the 
very wide range of value, that have been measured particularly at lower speeds, 
owever, s approaches a limiting value of less than .09, corresponding to 
tan s . 16 , for winds In e«cea. of 40 knots. Fro. Figure 22 It can be seen 
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Figure 25. Dependence of 
a* on wlndspeed at 
9.0 GHz. 

• Measured values 
— Theory 

(From Barrick, 1974) 



Figure 26. Dependence of o® on 
wlndspeed at 13.9 GHz. 

• Measured values 
•— Theory 

(Measurements from Jones 
et al, 1977) 
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Figure 27. Dependence of sea 
surface slope on wlndspeed. 
(From Barrick, 1974) 
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0 

(°) 

(db) 

Error in 
a° (+db) 

0 

(°) 

(db) 

Error In 
o° (±db) 


3ms“^ 

Hori 

zontal Polarization 

Ver 

tical Polar i 

zation 


UPWIND 

1.26 

11.87 

0.72 

1.25 

11.51 

0.72 

- 

DOWNWIND 

0.82 

12.49 

0.76 

0.77 

12.16 

0.76 


CROSSWIND 

0.69 

13.44 

0.72 

0.94 

12.86 

0.72 

i 

6.5ms~^ 

Horizontal Polarization 

Van 

deal Polarization 


UPWIND 

1.10 

11.19 

0.67 

0.89 

11.17 

0.67 


DOWNWIND 

1.06 

9.67 

0.70 

1.05 

9.39 

0.69 


CROSSWIND 

1.34 

11.14 

0.67 

1.37 

10.77 

0.64 


13.47ms“^ 

Horls 

sontal Polarization 

Vert 

deal Polarization 


UPWIND 

0.12 

9.12 

0.59 

0.12 

9.04 

0.59 


DOWNWIND 

0.12 

9.51 

0.59 

0.11 

9.41 

0.59 


CROSSWIND 

1.04 

9.26 

0.55 

1.06 

9.34 

0.55 


15ms~^ 

Horizontal Polarization 

Vert 

ical Polaris 

lation 


UPWIND 

1.33 

10.41 

0.75 

1.28 

9.96 

0.74 


CROSSWIND 

1.42 

9.69 

0.77 

1.76 

9.55 

0.74 


24ms“ ^ 

Horiz 

ontal Polarization 

Vert 

leal Polariz 

.at ion 


UPWIND 

0.41 

6.93 1 

0.54 

0.35 

6.67 

0.53 

I 

1 

I 

DOWNWIND 

1.03 

6.58 1 

0.54 

1.05 

6.54 

0.54 

( 

CROSSWIND 

0.30 

6.38 I 

0.53 

0.27 

6.38 

0.53 

I 

J 


6 is the mean angle of incidence in degrees. 

Table 8. Normalized Backscatter Cross-Section at 13.9 GHz (From Jones 
et al. , 1977). 
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that o®/|r( 0)1^ should vary from a minimum value of about 11 db up to perhaps 
25 db. The measurements of o® by Genest in Figure 25 and Jones et in 
Figure 26 demonstrate that such a minimum of is Indeed approached with a 
value close to that given by Barrick's theory. 

Meteorological data indicates that surface wind speeds in the range 3 
to 12 ms~‘ overwhelmingly predominate so that, in the absence of foam, at 

any frequency will generally vary by only a few db although occasionally it 
will substantially increase in exceptionally calm conditions. 

No measurements have been made to determine directly the effect of foam 
on sea surface backscatter cross-scatter cross section. However, radiometric 
observations of the microwave emisslvity, e, of sea surface show that at 9.4 
GHz (Williams, 1971) and at 19.3, 22.2 and 31.2 GHz (Nordberg et al . , 1971) 
e Increases from .4 to nearly 1 as the density of foam coverage increases. 

The distribution of bubble sizes as determined by Monahan and Zletlow (1969) • 
Figure 28, reaches a maximum of .5 mm radius and falls off rapidly so that few 
have radii in excess of 2 mm. Resonance effects can therefore be ruled out 
for the frequency range of interest here and e can be expected to be similar for 
all the pressure sounder channels. Since, e = 1 - |r( 0)P the presence of 
foam will reduce a°, at worst, in proportion to the fractional foam coverage. 
Ross and Cardone (1974) have related wind speed to foam coverage but the 
curve of Figure 29 smooths out the large fluctuationc- expected since the state 
of development of the sea also affects foam coverage. For wind speeds up to 
15 ms”^ there is usually less than 5% coverage but as much as 30% coverage has 
been observed with winds of 25 ms"^. Under these extreme conditions back- 
scatter cross section can thus be reduced by about 1.5 db and this could account 
for the deviation from theory of Jones 24 ms"^ measurement of in Figure 26. 

The experimental and theoretical evidence reviewed in 10.2 through 10.5 
substantiates the earlier assumption of an average value of at no‘ mal 
Incidence in excess of 5 db with a frequency dependence expressible by a 
second order polynomial. 

10.6 The Effect on Pressure Measurement of Higher Order Frequency Dependence 

in o® 

An Instrument with operating frequencies satisfying conditions 7.9 gives 





a measurement of pressure which is not affected by the second order frequency 
dependence of o“ given in Section 10.3. Since conclusive measurements of a° 
have not been made over a wide frequency range there remains the possibility 
of 3rd and Ath order dependencies of o° on frequency. In this case a° would 
have the form 

0 ° ■ exp (a + bv + cv^ + dv^ + ev**) 10. A 

For a particular set of operating frequencies It is possible to quantify the 
change in o^ which would give a change in the Index S, equivalent to a 
change in the surface pressure of the nominal design accuracy, 3 mb. 

This calculation has been carried out for the six frequency system of 
Table 2, line 4 and the measurement technique Is found to be relatively in- 
sensitive to high order frequency dependence. To give an error In S equiva- 
lent to a 3 mb change In pressure, o° must change by 1.5 db between 25 and 75 
GHz due to 3rd order dependence alone or by 0.75 db for the same frequency 
range due to Ath order dependence aloiie. 

10,7 o(6) Away from Normal Incidence 

TJhile an ability to remotely measure pressure along a subsatellite path 
is of great significance, it would be of even greater value if the two-di- 
mensional pressure field could be mapped by making off-nadir measurements. 

The limiting factor is the reduction of a (6) at angles away from nadir and 
this is most severe for calm sea conditions. Figures 22 and 27 allow a quan- 
titative assessment of these limits to be made. 

Surface wind speeds of a few ms ^ and their associated sea states fre- 
quently occur but only rarely does the wind speed fall below 1 ms ^ and the 
corresponding sea state occurs less often since the calm wind must prevail 
for many hours before the sea surface becomes calm. We can therefore expect 
B 0.04 frequently and s^ ■ 0.01 rarely. For these tan"*^ s is 10® and 5® 
respectively and in Figure 22 the corresponding values of o® are plotted as a 
function of angle of incidence. 

Off nadir measure will be limited by the reduction in o'’ and the conse- 
quent reduction in signal to noise which can be tolerated. Presently envisaged 




designs allow for a combined loss due to and attenuation by cloud and 
rain of about 10 db. Since calm sens are usually accompanied by little 
cloud, a minimum value of o ■ -5 db should still enable useful mcasuccmcnts 
to bo made. From Figure 22 It is apparent that measurement could be made 
at angles of up to 15** from the vertical except on occasions of n calm sea. 

For these low angles of Incidence the value of o® Is still dominated 
by tho gravity wave structure so that frequency variation due to Bragg re- 
sonance docs not occur (Gulnard et al. 1971) but frequency dependence through 
the Fresnel reflection coefficient applies as before. 

10.8 The Statistics of Sea Surface Reflectivity 

It Is common practice when considering the statistical variation of a 
microwave signal reflected by the ocean in a particular direction to assume 
that the surface is composed of many randomly— distributed, specularly re- 
flecting facets with appropriate orientations. The experimental results of 
Welssman and Johnson (1977) confirm that the distribution of wave heights is 
approximately Gaussian. Then, provided that the radar footprint Is larger 
than any coherence length in the sea-surface wave field, the return signal 
of a bearawidth-limited monochromatic radar is the (voltage) sum of a very 
large number of echoes, the in phase, x, and quadrature, y, components of 
which have Independent Gaussian distributions. 

The total power received is 

I(t) “ x’(t) + y^(t) 10.5 

where X ■ Ex and Y » Ey and the summations are over all the reflecting 
elements. The distribution of power then has the form 

which is the Rayleigh probability distribution for which the mean is I and 

o 

the r.ro.s. deviation 


is also Iq. The very large fluctuations In signal amplitude that this 
Implies Ifl known to observers as Rayleigh fading, 

It Is proposed to measure the ratio of the retutn fslgnnls at a pair 
of frequencies and If the Instantaneous reflectivities at the two frequen- 
cies are related a single measurement Is sufficient. However « Uelssman's 
dual-frequency radar experiment shows that coherence of the signals extends 
over a frequency separation of less than 40 MHz, Alt4)Ough the frequency 
coherence width Is greater for a space platform It does not approach the 
several GHz separation which la required for the pressure sounder. Conse- 
quently it Is necessary to effectively measure the average value of the 
normalized sea-surface backscatter cross section. 

To obtain an accurate value for the average reflected power (or 
equivalently, the average o^) which has the probability distribution of 
equation 10.6 a large number of Independence samples N must be averaged. The 
standard deviation of the average Is Io/>^ and the fractional accuracy is 
i/v*N. In Section 7.4 a fractional accuracy in the signal processing was 
Identified which was related to the design accuracy of the Instrument, 

6^3 mb/lOOOmb. Sufficient Independent samples must therefore be obtained 
to ensure that the contribution to the error in pressure measurement from 
the teflectivity statistics Is below this value. For a fractional accuracy 
of between .01 and .001 then N must be between 10** and 10®, The methods by 
which this number of independent samples can be collected are discussed in 
the following Section. 
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XX. Syatem Design and Measurement Errors 

The parameters svniXabXe for the design of a Microwave Pressure Sounder 
were identified In Section 3. The most Important are the operating fre- 
quencies and In Section 7 a detailed analysis was presented for the selection 
of the optimum frequencies. In this section we establish criteria for de- 
ciding values of the other parameters. The analysis is dominated by the 
accuracy requirement. The r.m.a. pressure error of + 3mb has been trans- 
formed using the pressure sensitivities of the chosen operating frequencies 
to a measurement accuracy specification, equation 7.22. The two major 
sources of error whic' interrelate all of the parameters are the bca-surface 
reflection statistics and the slgnal-to-noise ratio. It Is necessary to 
determine a set of system design parameter values so that the total error 
is within the specification. 

The size of the antenna will be determined largely by the need to 
collect many independent samples and this fixes the receiver bandwidth. 

The integration time and orbit height are also related to the antenna 
radius through the sea-surface error but are further limited by coverage 
and resolution requirements. The final parameter, transmitter power, has 
no optimum value but should have as high a value as possible to reduce 
the signal-to-noise error. 

The sea-surface statistical problem and its consequences for system 
design are considered first. 

11.1 Collection of Independent Samples and Signal Bandwidth 

The antenna size to be used can be determined from a consideration 
of the methods for collection of the independent samples required for an 
accurate measurement of o° as discussed in Section 10.8. Statistically 
Independent samples are produced in the following ways: 

1. By a change in the reflecting surface due to Its own motion. 

2. By illuminating a different reflecting surface due to motion 
of the satellite. 

3. By a (small) change in the radar frequency. 

4. By illuminating a different reflecting surface by swinging 
the antenna beam direction. 


M 


5. By using multiple antennae. 

Methods 1 and 2 occur naturally and we must examine whether they 
provide an adequate number of Independent samples within a reasonable inte- 
gration time. The other methods require some appropriate action and may be 
used to increase the number obtained naturally or to provide some advantage 
In system design which outweighs the disadvantage that these Imply for sys- 
tem complexity. The dctallf of each alterratlve are given In the following 
subsections. 

It Is convenient to Illustrate the various system options ns they are 
discussed in tliis section with numerical examples and to base these on a 
single set of parameter values. We shall use the following set and refer to 
it later as System A: 


Speed of satellite 

V 

7.2 km s"^ 

Antenna radius 

r » 

25 cm 

Orbit height 

h « 

500 km 

Integration time 

t « 

10 sec 

Duty cycle 

D “ 

.167 


The analysis is primarily intended to determine the antenna radius and the 
above value will only be used as a reference figure when this is not the 
object of the calculation. The orbit height and integration time given here 
are reasonable intermediate values and further discussion of coverage and 
resolution for a range of values will be given in Section 14. The duty 
cycle represents the fraction of the total time that each transmitter is 
switched on. We also chose, for illustration, an accuracy of 0.5% in each 
channel which, for the six frequency schemes of Table 2, represents a 1,5% 
accuracy overall or ±2mb rms error in pressure. 

11.1.1 Surface motion 

To change the reflecting surface sufficiently for statistical inde- 
pendence of the echo signals the linear changes in the positions of the 
specularly reflecting facets must be of the order of half the wavelength. 
The worst case conditions which limit reliance on this method for the pro- 
vision of an adequate rate of collection are for a calm sea where it is not 
unreasonable to expect speeds as low as 0.5 ms In these circumstances 
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tho coherence time of the 30 a surface for X/2 “ 6 mm (the longest wavelength 
presently envisaged) Is about 12 ms. Thus it would take -v 480 s to obtain 
the number of samples (N) needed for 0,5% measurement accuracy (1 //n) In a single 
channel and the integration time for tho six frequency instrument is totally 
unsatisfactory. 


11.1.2 Satellite Notion 

11.1.2.1 Spatial Coherence Length 

Satellite motion means that the reflecting surface is constantly 
being changed . The displacement of the satellite for statistical inde- 
pendence of two successively reflected signals may be obtained by application 
of the Van Clttert-Zernlke theorem (Born and Wolf, 1975). This gives an 
expression for the complex degree of coherence between signals at two points, 
Pi, P 2 , in a plane illuminated by an extended quasi-monochromatic source. For 
a uniform circular source of radius R there is complete Incoherence of the 
signals when the separation of Pi and P 2 is 

d = 0,61 hX/R 11.1 

where h is the distance between the source and the plane of P 1 P 2 . Applying 
. wSnt situation, the illuminated area of sea surface is the 
lich approximates well to a continuous distribution since 
^ large number of reflecting facets. The radius of this 

circulr.r uoarc 3 related to the angular beam width m of the radar antenna 
by 

to = R/h 11.2 

The angular beam-width may be approximated by the angle to the first minimum 
of a diffraction limited aperture (the antenna of radius r) 

to ■ .61 X/r 11.3 

Combining these three equations gives d = r. However we can expect 
the phase changes between Pi and P 2 to be double that implied here since the 
phases in the extended source are affected by the motion of the radiating 
antenna. So that the separation of points for complete Incoherence is 


d = r/2 


11.4 


Independent samples arc thus obtained by displacement oC the antenna by 
this jpatial coherence length of half the antenna radius. 

The limitations which make this an approximation to within a small 
numerical factor are; firstly that the extended source is not uniform but 
has a Bessel function distribution and the radius used In equation 11.2 is 
that of its first minimum; secondly the signal is detected not at a point 
but is the sum over the antenna area. A more rigorous analysis of the 
problem by Peckham (1975) takes these factors into account and gives the 
spatial coherence length as 0.89 r for a circular antenna. 

With this value of the spatial coherence length, 1.5% accuracy is 
obtained in a time 0.89 r / [VD(0.5 x 10 * 2 ) 2 ] seconds , which for System A 
is 7.7 s. Thus the rate of collection of samples by the natural satellite 
displacement overwhelms that due to changes in the sea surface with time. 
When other methods for increasing the collection rate are not used, the 
antenna radius is fully determined by the required fractional accuracy, F^, 
and is given by the relationship 

t V D Fa2 11.5 

^ ~ 0.89 

It is apparent from the above numerical example that for the six frequency 
sets of Table 2, is large enough for adequate collection rates to be 
achieved by this natural process with a reasonable antenna size. 

Of course, it is not essential to use a circular antenna. The most 
favorable arrangement would be to have a small antenna dimension alottg the 
direction of satellite motion, which ensures a high sample collection race, 
and a large transverse dimension to maintain a large receiver aperture and 
hence small slgnal-to-noise errors. However, the antenna dimensions deter- 
mine the surface spot size so that ground resolution considerations will 
also influence the ultimate choice of aperture shape and size. 

For a non-circular antenna the spatial coherence length must be re- 
examined. Peckh.m (1975) shows that a rectangular aperture of side 2a in 
the direction of motion has a coherence length of 1.08a and it can be ex- 
pected *-hat for an elliptical antenna the coherence length is between this 


«* 


and the cLreular valtie Ic about o(|ua 1 t o tlu‘ .si*mi ^in i nor axis. 

» 

11.1.2.2 Band wid til 

A Microwave Pressure Sounder relyluR on satellite motion to provide 
independent samples would be operated in an intermittent c.w. mode. The 
return signal is ^ 130db below the transmitted power and the necessary 
isolation is not possible to allow simultaneous transmit and receive operation. 
Consequently a fixed frequency signal would be transmitted for the few 
milliseconds required for the sigucti to return after reflection at the ocean 
surface. The return echo has a larger bandwidth than the transmitted 
signal because of Doppler frequency shifts produced by components of the 
satellite velocity over the angular spread of the beam. The velocity com- 
ponent along the 3db beam direction is ±0.25 VX/r so that the Doppler shift 
is ±0.5 V/r and a total 3db receiver bandwidth 


Af = V/r 11.6 

is needed to match the natural bandwidth of the return signal. For System A 
this has a value close to 30 kHz. 

Satellite pitch instabilities will cause deviations of the antenna 
boreslght from the nadir direction and will consequently offset the band of 
received frequencies from its symmetrical arrangement about the transmitted 
frequency. The specified inertial pointing accuracy of Spacelab is ±0.5° 
and since this is of similar magnitude to the beamwldth the frequency offset 
would be comparable to the bandwidtli. Methods are available for impvwing ' 
the pointing accuracy of the platform and hence reducing the offset if this 
is desirable. However, this Is not a restricting problem so long as the 

overall pressure sounder system Is designed to accommodate such bandwidth 
offset effects. 

11.1.3 Frequency Sweeping 

A combined consideration of the sea surface statistical error, the 
slgnal-Lo-nolsc error and surface resolution suggests that significant im- 
provement in accuracy Is possible if the natural rate of collection of in- 
dependent samples can be Increased by one of the techniques (3-5) identified 
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above. Detailed consideration must therefore be given to the potctinl a. - 
vantages and engineering fenslbllitv of each alternative so that it ma, be 
decided which if any should he implemented. We begin in this subsect 
with the frequency sweuplns technique. 

The method of Section 7 for choosing the operating frequencies 
has assumed that each is monochromatic. Small deviations from these values 
can be tolerated without destroying the designed insensitivity o t e com 
bined signal to unwanted background, temperature and water vapor 
However, if a sweep of 'v- 100 MHz is required then the effects of t ese a mo- 
spheric parameters on the accuracy of the pressure measurement must be re- 

examined . 

The fundamental reason why sweeping over a range of frequencies cap 
provide a number of independent samples is that signals at a pair of dif- 
ferent, but close, frequencies will have nncorrelated received echo ampli- 
tudes after reflection at the sea surface if the frequency separation is large 
enough. To increase the sample collection rate over that ptodne^ by satel- 
lite motion the sweeping must be accomplished within the time taken for 
satellite to be displaced by the spatial coherence length, le within e time 
t. - r/V ( - 34 PS for System A) . The following analysis detetmlnes how much 
Che frequency must be changed to produce a statistically independent sample. 

11.1.3.1 Frequency Coherence Width 

Weissman (1973) derives a correlation function R(Ak) for two (voltage) 
signals reflected from the ocean at normal incidence when the frequency se- 
paration of the two radar signals is Av = Ak, (c is the velocity o 
electromagnetic waves). The correlation function has a magnitude 


lR(Ak)l * exp [-2o^(Ak)q 


11.7 


|(Ak)he32_| = 
^ '*'1 5 .544 i 


1 . 

2 


where o is the r.m.s. ocean w.rve height .ind the other symbols have been pte- 
vlousl, defined. The two factors which contribute to the cotrelntlon function 
derive independently from the wnve height statistics and the antenna pattern 
function. This analytic expression for R(Ak) was obtained by Weissman by 




Integrations over the wave height distribution and the antenna bcamahape. 

The Integrations were facilitated by assuming a Gaussian distribution of 
wavehelghts and a Gaussian beamshape. 


We require the separation, w^ ® Wj^, of two frequencies such that 
the two received signal powers are statistically independent. This we define 
as the frequency coherence width and it may be expressed in terms of Welss- 
man's correlation function by 


w 


k 


lR(Ak )|2 d(Ak) 


11.8 


Approaching the problem in a slightly different way Peckham (1975) 
calculates Wj^ for that part of the correlation function due to the beam- 
shape. He uses a Bessel function antenna pattern function and by numerical 
integration obtains 


w, = 0.5672 irk^r^ 11.9 

k -~jp- 

where k = 2 ttv/c. 

Alternatively, Wj^ may be obtained from the beamshape factor of 
equation 11.7 thus 

10 


Now, the 3db beamwldth, 63 , for an antenna beam shape given by G(6) «2Ji (kr9) ^ 
is 3.2326/kr, so that 

w, • .5305 irk^r^ 11.11 

h 

and this is in good agreement with equation 11.9 above. It is readily shown 
by substituting the coherence width, equation 11 . 10 , into equation 11.7 that 
for a frequency separation of Wj^ the power correlation coefficient ] R(Ak) p 
has the value (1 + = 0.092 when the ocean is smooth (o << 1 ). 


w. 


/ 


d(Ak) 


I 5.544 I 


= 5.544ir 


h0,2 


11 
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11.1.3.2 P Surface Effect on Frequency Coherence Width 


The effects of sea surface statistics on the frequency coherence width 
may be determined by repeating the integration, equation 11.8 with both 
factors of equation 11.6 Included in |R(Ak)p. The new coherence width 
is best expressed in terms of the previous Wj^, equation 11.11. With the 
help of standard Integrals we obtain 


and the function 


1 

exp 


\ f 


.'IT ' , 


11.12 





11.13 


o 

can be obtained from tables of the error function. To facilitate interpreta- 
tion of this result, equations 11.11 and 11.12 may be written explicitly as 
frequency coherence widths, Wy and w^, in terms of the frequency v: 


and 


w = 10.5 


V ch 


i 

1 

7 > 

> 

II 

:> 

f / 4ow,, \ 

exp[(— v) J 


11.14 


11.15 


wl 

Then values of the ratio of frequency coherence width are tabulated for 

A \ 

a range of values of the argument ( ie for a range '^of sea states) in 

Table 9. 


4ow^ 

c 

.1 

.2 

.5 

1 

2 

5 

10 

20 


.99 

.80 

.62 

.43 

.25 

.11 

.056 

.028 


Table 9. Effect of Ocean Surface Statistics on the Frequency 
Coherence Width (see text for explanation). 
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The tabic shows that roughness of the sea surface reduces the fre 
quency coherence width, le wj, 1« always less than w^,. A -smooth- sea, 
where 4ow /e i Is sufficteiU deflnlllmi of a smooth ;.ea lor pr.'sent pur- 
poses, ims little effect on the coherence width. However when the argument. 
4ow /c, is large the ratio, w^w^^ approaches tlie value ( 

the frequency coherence width is determined principally by the sta surface 
effect rather than the. antenna beam shape. 

For a given set of experimental parameters the relative importance of 
surface height statistics and beam shape on the frequency decorrelation can 
now be determined. For System A above with v » 50 GHz the frequency coherence 
width w = 10.9 MHz and if o - Im (a moderately rough sea) (40 w^/c) = 0.14. 
Thus the beam shape is the dominating factor. The 14 GHz correlation radar 
used by Weissman and Johnson (197?) to determine wave heights had a Im dia- 
meter antenna and was operated from an aircraft at 3 km altitude. In this 
case w^ = 4.8 GHz but since (40 w,/c) = 64. at o = Im. the frequency co- 
herence width, when the sea surface effect is included, w^ = 42 MHz. As ex- 
pected the surface effects dominate. The main reason for these two example 
systems being at either end of the scale is the difference in altitude bet- 
ween the satellite an- -aircraft platforms. 


In a swept frequency system, suitable for pressure sounding from a 
satellite and therefore, as presently envisaged, with parameter values not 
very different from System A, the beam shape will be the important factor 
and equation 11.14 is a satisfactory approximation for the trequency co- 

herence width. 


11.1.3.3 Implementation and Bandwidth 

A frequency swept system requires a larger value of bandwidth than 
the fixed frequency scheme. The range of frequencies being received at any 
instant is determined by the rate at which the frequency is swept and the 
time delay between the signals reflected from immediateJv below the satellite 
and from the edges of the antenna beam. If the sample collection rate, is lo 
be increased by a factor m, and thus the sea surface statistical error re- 
duced by m"'^ then the total sweep 


Av. 


= mw. 


Hz 


* 


11.16 
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must be completed in a time r/V and the sweep rate is 

6v„ • m 10.5 v^rV Hz s“^ 11.17 

“ch 

A reasonable minimum value of m is 9 to reduce the statistical error by 1/3. 

This gives for System A 

Av - 100 MHz and 6Vg = 3 MHz ys"^ . 

These values are rather large and would necessitate development of suitable 
millineter wave sources with either current or varactor tuning. In addition 
the large total sweep means that the frequency selection procedure should be 
revised to ensure retention of the insensitivities previously specified. 

It can be readily shown that the minimum bandwidth required is 

Af ® A.2 m V 11.18 

r 

where V/r is just the bandwidth for a fixed frequency scheme. The factor of 
4.2 in this equation makes the signal-to-noise for any swept system at least 
3db worse than for a fixed frequency system. Specifically, when m “ 9 the 
bandwidth is 1 MHz and the S;N is 8db worse. 

These values assume that the frequency swept mode is used to full 
advantage by having a tracking local oscillator such that the IF band and 
the received band of frequencies are the same. For these to coincide within 
1% of the bandwidth the implied phase sychronization of the local oscillator 
sweep modulation is equivalent to a knowledge of the satellite altitude within 
(c/24 Wv) meters Imfor System A). Continuous operation over the ocean, 
where the altitude changes only slowly, does not present a severe sychroni- 
zation problem and can probably be accomplished with appropriate processing 
of the pressure sounder signals. However, after each land-to-ocean coastal 
crossing the local oscillator phase will need to be reiocked. 

The use of a tracking local oscillator with a larger receiver IF than 
the minimum, 11.18, does not significantly reduce the synchronization problem. 
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A simpler mode of operation is to use a fixed frequency local oscillator and 
a receiver bandwidth of Av„ but this will reduce S:N even more and is there- 
fore probably not viable. 

Besides providing, an increased sample collection rate, frequency 
sweep Inp, also gives a considerable degree of range resolution. However, this 
does not constitute a valid reason for choosing a swept frequency system 
since backscatter from aerosols is negligible when cloud attenuation is small 
enough to allow operation of the pressure sounder. 


As alternatives to continuous sweeping of the frequency, step changes 
equal to the frequency coherence width may be made to obtain the independent 
samples or nanosecond pulses may be used. These techniques do not appear to 
offer any additional advantages and are no easier to implement. 


11.1.4 Antenna Scanning 

An independent sample can be obtained with a fixed frequency radar by 
swinging the antenna beam through half the angular beamwidth so that more 
than half the area of sea illuminated is changed. The scanning can be 
either along the sub-satellite path or transverse to it. The rate of col- 
lection of independent samples can be increased by a factor m over the 
natural rate by sweeping through m/2 beamwidths in a time r/V. Since r/V 33 ps 
for System A this implies scanning rates of several kHz. Mechanical scanning 
of the primary antenna at this rate is not possible and the problem is only 
slightly reduced if scanning with the feed or a sub-reflector is considered. 
Because six frequencies are needed, electronic scanning with a phased array 
is complicated unless several antenna are used and if a multiple antenna 
system is anticipated then it is probably better used as Indicated in Section 

11.1.5 below. 

Antenna scanning does not appear to be feasible at this stage. 
Nevertheless it is worthwhile to Identify some of the advantages of the 

technique. 

With antenna scanning a larger antenna can be used to advantage. For 
a swept frequency scheme a larger antenna ..lakes implementation more difficult 


by lucroaainn both thc> total p.weep raiipo and the swoop rate to maintain the 
rate of collection of indepondenl; samplon. In contrast, with a scanned 
antenna system, n larnt^r antenna requires to be scanned over the same total 
an>^le hut at a reduced sweep rate to maintain tlic collection rate and the 
teciinlques boeomes more feasible, A larger antenna Is desirable because of 
the Improvement In received signal strength and the reduction In noise 
through a narrower Doppler bandwidth. The consequent Improvement in S.N 
allows more accurate measurements under worse atmospheric conditions. 

The scanned ,antenna system is more complicated than the frequency 
swept scheme If it is to be used to full advantage because synchronization 
of the transmit and receive antenna scans must be used in addition to the 
control of the local oscillator frequency (for along track scanning). When 
used in this way there is also a degree of height discrimination. The 
tight tracking constraints can be relaxed with a reduction in the advantages 
by using antennae of unequal size and widening the receiver bandwidth. This 
leads to numerous alternative methods for implementing an antenna scanned 
system each of which has different consequences for the rate of collection of 
independent samples, the S:N ratio and the height discrimination while each 
Imposes different constraints on the system design in terms of sweep rate, 
precision required in sychronized tracking and local oscillator frequency 
control. All of these quantities can be readily determined for a specific 

design. 


11.1.5 Multiple Antennae 

The sample collection rate can be increased by a factor m by using 
m multifrequency antennae, each with its own IF chain and square law de- 
tector. Although this has the advantage of redundancy and can be used in 
combination with any of the previous techniques, size, weight and cost con- 
siderations suggest that the concept is only of value when the antennae are 
small and an increase in the collection rate of less than ten is desired. 

11.2 Signal to Noise Errors 

The sources of noise which contribute to the error in pressure mea- 
surement and must therefore be taken into account are: 


f 
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of null-balancing techniques should be utilized for control purposes. Al- 
though a comparison of the received and transmitted powers could be made 
at each frequency It is probably more satisfactory to control the ratio of 
transmitted powers at a pair of frequencies, P,j(vi)/P,j,(vp) , so as to 
equalize the received echo power in the absence of atmospheric attenuation. 

Fluctuations In the gain of the IF amplifier can be overcome by 
using the same IF chain and detection system to amplify and ratio the i c- 
celved signals of a frequency pair . 

Although a A— frequency, 3-ratio scheme Hiuy appear simpler than a 
6-frequency instrument, the system architecture detailed above means that 
the power outputs at all four frequencies must be related and the detector 
must accept all four return signals. Consequently a 6-frequency scheme with 
three separate ratlometers is an easier option to implement. 
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12.1 Fixed-Frequency Microwave Pressure Sounder 
12.1.1 Primary System Design Parameters 

The operating frequencies and other design parameters of Section 3 

must be supplemented by the system efficiencies, n^ and Dq, the satellite 

i K 

speed, V, and the duty cycle, D, to complete the Instrument specification. 
Table 10 gives the values of these basic parameters for a fixed frequency 
system operating at two altitudes. Since a rectangular antenna is speci- 
fied, the radius has been replaced by the overall dimensions. 


V. 

1 

Operating Frequencies 
(GHz) 

29.2555 A4.80 67.51 

36.5555 52.80 73.01 

a, 6 

Indices 

-1.60, 1.00 

Pt 

Transmitter power 

2W 

H9B9ili 

Feed efficiencies 

0.85 

(r) 

Antenna dimensions 

20 X 150 cm 

t 

Integration time 

12 s 

D 

Duty cycle 

0.166 

h 

Satellite altitude 

500 

800 km 

V 

Satellite speed 

7.6L 

bibhedhhhii 

Af 

Receiver bandwidth 

76.1 

Ih . 5 kHz 

2h/c 

Pulse length 

3.33 

5,33 ms 


Table 10. rixed Freepumey System ])esi; n Parameters 


















12.1.2 Derived Performance Data 


The parameter valuea of Table 10 enable the instrument performance 
to be evaluated. Those data which are not dependent on atmospheric surface 
conditions are given here in Table 11. 


Sensitivity of Index S 
to pressure change 

1. 


0.74% per mb 

Sensitivity of pressure 
measurement to: 

a) Background variation of the 
form (a + bv + cv^) 

b) Surface temperature 

c) Total water vapor 

From 

Frequencies 

negligible 

800°C per mb 
8 gm cm“^ per mb 

R.m.s. atmospheric 


From numerical 

± 0.4 mb 

variability error 


simulation 


Signal processing accuracy 

From 1 above 

0.73% 

required in each channel for 



I ± 3mb r*m.s. error 




Number of 

independent samoles 

N = 

1.4 X 105 

1 per channel in time t 


1.08r 


Accuracy per channel 


0.27% 

Total statistical error 

Sum over all 

± 1.1 mb 




channels 


Beamwidths and Surface Resolution 

1 

Channel ' 

Frequency 

Wavelength 

3db 

Ground Resolution (km) 


(GHz) 

(mm) 

Beamwidths (°) 

h“500 km h»800 km 

1 

29.2555 

10.25 

1.30 X 0.17 

136 X 6.8 162 X 10.9 

2 

36.5555 

8.21 

1.04 X 0.14 

126 X 5.5 145 X 8.7 

3 

44.80 

6.70 

0.85 X 0.11 

118 X 4.5 133 X 7.1 

4 

52.80 

5.68 

0.72 X 0.096 

113 X 3.8 125 X 6.0 

5 

67.51 

4.44 

0.56 X 0.075 

107 X 3.0 115 X 4.7 

6 

73.01 

4.11 

0.52 X 0.070 

105 X 2.7 112 X 4.4 

Effective along-track ground 



resolution, all channels 


85 km 80 km 


Table 11. Fixed Frequency System Performance 
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The signal processing accur'iru 
tac.ion I'^dcToTo.; 

tion. " required from the instrumentn- 


cnannel dependent ground resolution in this table t m 
-rtac. are., ,„el„aa<l „Uhi„ the first sere of th ’ 

---- - oT:::;::::'"" 

pathln'th "Thri?thtefl?t'riff! u "T 

and Is the same tor all channels. ^ 

12.1.3 Total Error Calculations 

The values of two components of the total error t-u 
variability error E emi ^ atmospheric 

-ave heen liven I: r^h^^e l ' ^SS- 

third source of error si '1 »>’atem, the 

pandent on atmospheri: and ocel°:u";::et„d^:lte:“t:t\:l!:\" 

!Z"e^"'' d”" S'- a<!“atlons 11.19 endTrir"'' 

can be used to evaluate this error P , k . ’ 

quencv la from enuntlon 10.1 


Atmospheric I 
1 Transmission Im 

Channel 

ceiver noise figure 
Pg = 1013mb 
cra~2 

^ , Cloud 1 grom”"^ 

1 

7.0 

-0.18 

-0.47 

-1.29 

2 

7.0 

-0.33 

-0.52 

-1.99 

3 

7.5 
-0.88 
-0.72 
-2.97 j 

4 

8.0 

-9.14 

-0.99 

-4.09 

5 

8.5 

-7.31 

-1.59 

-6.60 

6 

8.5 

-1.68 

-1.86 

-7.69 

Backscatter 
Cross-Sec tior 

0 Very calm sea 
Moderate sea 
Very rough sea 

0(14'’) Very calm sea 
Moderate sea 
Very rough sea 

' 15.48 
10.48 

6.48 

-10.52 

2.48 

4.48 

15.28 

10.28 
6.28 

-10.72 

2.28 

4.28 

15.07 

10.07 

6.07 

-10.93 

2.07 

4.07 

14.87 

9.87 

5.87 

-11.1 3 

1.87 

3.87 



14.51 

9.51 

5.51 

-11.49 
1 .51 

3.51 

1 :4.38 

9.38 

5.38 

-11.62 

1 . 38 

3.38 


Table 1:’. 

TypicnJ 

Atnio,‘ipherlc and 

Surfai‘0 

Data (d 

h) 


In Table 13 the total measurement error in mb Is shown for both 500 km 
and 800 km orbit altitudes and for various atmospheric and sea surface states. 
The total Is obtained by adding the three components Egg, E^y and Eg^j as in- 
dependent random errors. It can be seen that the r.m.s. error with this in- 
strument is 1.2 mb for the most frequently occurring conditions but, as 
expected, it will not be able to provide sufficiently accurate results in 
heavy cloud or away from nadir over very calm sea. The 500 km orbit is 
slightly better than the 800 km orbit but the latter may be preferred 
because in this case measurements 14° from the nadir are 200 km from the 
sub satellite path compared with 125 km away for the 500 km orbit. Curvature 
of the Earth has been neglected, thus 14° from nadir is equated with 14° from 
vertical - an error of 1 . 8 °. 



Atmospheric Conditions 

Sea State 

Surface 

Pressure 

1013mb 

1013mb 

1013mb 

990mb 


Water 

Vapor 

2.6 gmcm”^ 

5.2 gra cm“^ 

3 . 5 gm cm“^ 

3.5 gm cm"^ 

500 km 

Uniform 

None 

0.5 km 

1 km 

2 km 

Orbit 

Cloud 


0.5 gm m”^ 

1 gm m”^ 

1.5 gm m“^ 

Very calm 

Nadir 

■M 

1.17 

1.17 


14° 


3.23 

>5 


Moderate 

Nadir 

1.17 

1.17 

1.17 

1.39 


140 

1.17 

1.18 

1.21 

>5 

Very rough 

Nadir 

1.17 

1.17 

1.18 



140 

1.17 

1.17 

1.18 


800 km 
Orbit 

Very calm 

Nadir 

— 

1.17 

1.17 

1.32 

14° 


>5 

>10 

>25 

Moderate 

Nadir 

1.17 

1.17 

1.18 

2.23 


14° 

1.19 

1.23 

1.42 

>10 

Very rough 

Nadir 

1.17 

1.18 

1.20 

4.92 

14° 

1.18 

1.19 

1.28 

>5 


Table 13. Total r.m.s. Errors (mb) With Fixed Frequency MPS 
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12,1.4 Hardware Iroplemontation 

The system specifications of 12.1,1 have been used to define an out- 
line hardware Implementation of the Microwave Pressure Sounder. The in- 
strument comprises three ratlometers which are of similar designs but differ 
in their operating frequencies. The block diagram of a two-frequency ratio- 
meter is shown in Figure 30. Each ratiometer has the following subsystems: 
power supply and frequency control; transmitter; RF output power equalization; 
receiver and ratiolng electronics. A single multifrequency dual-antenna sub- 
system serves all three ratlometers while common electronics Include timing 
and modulation control, signal processing and satellite power and data inter- 
faces. A pointing subsystem is also required. 

The instrumentation is based, with one exception, on existing techno- 
logy. The Gunn effect or Impatt sources used as the first stage oscillators 
of the transmitter subsystem have their frequencies crystal controlled. 

These are used to drive the injection-locked millimeter wave power oscillators. 
The devices at present available for this second stage do not provide suffi- 
cient RF power. However, suitable multi-diode Impatt devices with a gain of 
up to lOdb and an output power of 2W cw are currently being developed 
commercially. 

The two frequencies are transmitted alternately and this enables the 
output power equalization loop to utilize phase sensitive detector control of 
a ferrite modulator. The loop must provide equal output energies at the two 
frequencies over the Integration period or any sub-division of this for which 
the ratio is computed. Set level attenuators are used to accommodate asym- 
metries between the two branches and differential RF detector response. Ex- 
periments have shown that such a loop provides more than adequate equalization 
of the output powers. 

The most suitable antenna has a multi-frequency feed, similar perhaps 
to that used for the SMMR on SEASAT A. The corrugated horn feed on this 
Instrument has a high efficiency over a 6 to 1 frequency range. S' parnte 
transmit and receive antennae are preferred over a single antenna with b 
circulators because of the smaller loss. 



Figure 30. Two Frequency Ratiometer Schematic 
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It is possible to use a single mixer in the receiver but the large 
difference between a frequency pair means that its response to each fre- 
quency will be different. However, it is easier to match the response of 
separate mixer-preamplifiers than to do this with a single device. Con- 
sequently this arrangement is chosen for the front-end while the main IF 
amplifier chain and the detector are common. Automatic gain control with a 
long time constant is used to reduce the dynamic range required from the 
square law detector and the final ratioing is accomplished digitally. It 
has been demonstrated that the stability and accuracy specifications of the 
detection and-ratlometer electronics can be met. 

The modulation sequence for each of the separate ratiometers must 
allow for measurements to be made of the signal and noise at each frequency 
when the instrument is not transmitting. Operation of the transmitters 
is then in an intermittent c.w. mode with pulse lengths equivalent to the time 
taken for a signal to be reflected at the surface and returned to the re- 
ceiver, ie 2h/c seconds. A full cycle of operations takes six of these 
periods, two for transmission and 4 for detection. The large differences 
between the six frequencies will allow the three ratiometers to be operated 
simultaneously with synchronized transmit and receive sequences. Figure 31 
shows the proposed oscillator and detector modulation sequence. Dynamic con- 
trol of the repetition rate can be used to accommodate variations in the 
orbit height. It may be possible that, with sufficient filtering, noise 
measurement could be made at one frequency of a pair while the other was 
being transmitted. This would improve the duty cycle and reduce the surface 
statistical error by a factor of 0.82. 

Further specification of component performance, stability and toler- 
ances is required. Although the hardware design is based on available 
technology the RF power stability requirements for this instrument are not 
generally needed in radar systems and have therefore not previously been 
demonstrated. However, the use of phase sensitive detector control of the 
output powers and matched mixer-detectors should enable the specifications to 
be met. Some parts of the system have been investigated but further testing, 
particularly of the receiver subsystem, is needed. Careful preflight cali- 
bration and stability testa of the completed instrument will also be essential. 


« 



Figure 31. Modulation sequence of the two-frequency ratiometer. 
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2.2 Swept Frequency Microwave Pressure Sounder 
12.2.1 Primary System Design Parameters 

The exMple ot e swept frequency system used here for comperleon has 
e majority of Its design parameters the same as for the fined frequency 
example above. To these parameters ve must add the factor, m, by which 
frequency sweeping increase, the sample collection rate. This then defines 
the total sweep and the sweep rate, which are height and frequency dependent 
and the minimum receiver bandwidth. We choose a larger antenna to partially’ 
compensate for the worsening of slgnal-to-nolse by an Increased bandwidth. 
Heaver. If the total sweep Is to be kept within reasonable bounds « 100 MBs) 
t e Increase in antenna else must be limited. The complete system specif 1- 
ation is given in Table 14. 




(r) 


Operating Frequencies 
(GHz) 

Indices 

Transmitter power 
Feed efficiencies 
Antenna dimensions 


Integration time 
Duty cycle 
Factor 

Satellite altitude 
Satellite speed 
Receiver bandwidth 
Pulse length 

repetition time 


29.2555 
36.5555 
-1.60, 1.00 

44.80 67.51 

52.80 73.01 

2W ^ ■ — 

0.85 

28 X 150 cm 

12s 


0.166 


12 


500 

800 km 

7.61 

7.45 kms'^ 

76.1 

74.5 kHz 

3.33 

5.33 ms 

18.4 

18.8 PS ' 


Channel 

A. Total Sweep j 
Sweep Rate 

B. Total 


Sweep 8.81 
Sweep Rate .469 


22.0 

1.20 

13.7 

.732 


33.0 

1.80 

20.6 

1.10 


45.9 

2.49 

28.7 

1.53 


75.0 

4.08 

46.9 

2.50 


54.8 

2.92 


Table 14. 


Swept Frequency System Design Parameters 


12.2.2 Derived Performance Data 

The main difference in the performance of this system when compared 
with the details of Table 11 for the fixed frequency system is the improved 
rate of collection of independent samples provided by the frequency sweeping. 
The factor m and the antenna size both affect this and make the Instrumental 
-error attributable to the sea surface statistics ± 0.38 mb for this Imple- 
mentation. 

Detailed calculations have not been performed to find the sensitivity 
of the measurement to changes In the various atmospheric factors. However 
It Is reasonable to expect that this or a similar set of operating fre- 
quencies will result In a r.m.s. atmospheric variability error close to 
± 0.4 mb as found for the fixed frequency scheme. 

The effective along-track ground resolution, being the product of 
satellite ground speed and Integration time. Is not changed by the Increase 
In antenna size although this Increase Is along the direction of satellite 
motion. However, the reductions In beamwldths and smeared spot sizes mean 
that a greater proportion of the Integrated energy comes from within the 
effective resolution cell area. 

An abbreviated list of performance data Is given In Table 15. 


Rms atmospheric variability error 

± 0.4 mb 

Bms statistical error 

± 0.38 mb 

Effective along-track ground resolution; 
for satellite altitude: 

85 km 
500 km 

80 km 
800 km 


Table 15. Abbreviated Swept Frequency System Performance 
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12.2.3 Total Error Calculations 

The use of frequency sweeping sUows th- reduction In ses-surfece 
stetlstlcsl error, duo to the sweeping to be tredod for on Incrensc In the 
slgnsl-to-nolso error, due to the wider receiver bandwidth needed. The 
following result, show the effects of this trade-off for the exsnple systen 
specified above. It has been essuned that fullest sdvsnteges Is tshen of 
the swept frequency systes, by using the mlnlssim receiver bandwidth. 


The data of Table 12 has been used to calcuUte the slgnal-to-nol.e 
error. E,„. for the set of atnosperlc and surface conditions In Table 13. 

1. then combined with E^„ and Ejj to produce the total r.n.s. errors n 
si^surad pressure. The results presented In Table 16 are directly comparable 
with the values for the fixed frequency scheme In Table 13. 




Atmosp 

heric Conditions 



Surface 

1013mb 

— 

1013mb 

1013mb 

990mb 


Water 

Vapor 

2.6 gm cm“^ 

5 . 2 gm cm“^ 

3 . 5 gm cm"2 

3.5 gm cm”^ 

500 km 
Orbit 

Uniform 

Cloud 


0.5 km 
0.5 gm m“^ 


2 km 

1.5 gm m"^ 

Very calm 
Moderate 
Very rough 

Madlr 

1A° 

Nadir 

14° 

Nadir 

14° 

0.55 

>5 

0.55 

0.70 

0.57 

0.61 

0.55 

>10 

0.56 

0.94 

0.61 

0.73 


0.56 

>25 

3.84 

>25 

>5 

>10 

800 km 
Orbit 





Very calm 
Moderate 
Very rough 

Nadir 

14° 

Nadir 

14° 

Nadir 

14° 

0.55 

>10 

0.57 

1.22 

0.67 

0.88 

0.56 

>25 

0.61 

2.01 

0,87 

1.34 

0.58 

>25 

0.76 

4.14 

1.46 

2.65 

3.09 

>25 

>5 

>25 

>10 

>25 


Table 16. Total t.m.s. Errors (mb) With Frequency Swept MPS 
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Under favorable conditions the error has been halved to ± 0.6 mb 
but for some frequently occurring cloud and sea states the error is 
larger and the range of conditions for which the frequency swept system 
operates with sufficient accuracy is more restricted. It should also be 
noted that the limitations are more severe for the preferred 800 km orbit 
The improved accuracy has been bought at the expense of reduced coverage 
and increased system complexity. 

12.2.4 Hardware Implementation 

The frequency swept system will consist of three ratiometers each 
of which has the same subsystems as described in Section 12.1.4 for the 
fixed frequency scheme. The block diagram of one of these ratiometers, 
Figure 30, is also applicable although some of the components are used 
differently. 


Sweeping of the transmitter frequency must be accomplished with 
either current or varactor control. Extensive development will be neces- 
sary to produce sources with sufficient power and the required sweep 
characteristics. The local oscillator is to be swept over the same fre- 
quency range and accurate synchronization of its cycle must be. provided by 
altimeter type measurements which can be obtained from the pressure sounder 
signals. The modulation sequence is as in Figure 31 but each transmission 
period includes many sweeps through the range of frequencies. 


Successful operation of this frequency swept Instrument will depend 
upon careful design and calibration for the 20 to 100 MHz bandwidth used 
at each frequency. 
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12.3 Payload Physical Characteristics 

The additional data (approximate values) in Table 17 of size, weight 
etc. complete the payload spec if icat ions for either of the implementation 
olternatlves described above. 


Overall size 

1.5 X 0.6 X 0.5 m 

Weight 

< 100 kg 

Field of view 

5° X 0.5° 

Pointing angle 

Nadir; ± 14° (cross-track) 

Pointing accuracy 

± 0.2° 

Power 

<100 W 

Data Rate 

1 kbps 


Table 17. M.P.S. Physical Characteristics 


12.4 Platform Requirements 

The most suitable orbit for use of the Microwave Pressure Sounder 
as part of an operational veteorological data collection network is a near- 
polar, sun synchronous orbit at about 800 km. Potential missions on which 
to demonstrate the long-term capability of the instrument, prior to its 
operational application, are Oceansat, Cllmsat or the proposed Microwave 
Observatory. Shuttle/Spacelab sortie missions could be used for proof-of- 
concept flights with essentially the same instrumentation. 

Aircraft testing of the technique is possible only for a fixed fre- 
quency system. The difference in altitude between an aircraft and a satel- 
lite means that the frequency coherence width is too large to allow sweeping 
to be used to Increase the sample collection rate with the lower platform. 
Even with the fixed frequency system extensive modifications to the inte- 
gration time, bandwidth and modulation timing are required for aircraft 
operation. In addition, because the instrument is designed to measure the 
integrated effect of the total atmosphere, accurate measurements of pressure 
at the aircraft altitude would be necessary and a different set of operating 
frequencies might be required. Careful consideration would also need to be 



given to the atablUty requirements of the aircraft platform for meaningful 
meaBurements. Nevertheless, development of an aircraft instrument may be 
of value for its potential application to severe storm monitoring on flight 

paths through a storm. 
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13 , C ovfcTragR and Ref^.olutlo n 

The Mlcrownve PrcHSure Sounclnr ItiHtrument spec If lent lonn of Section 
12 enn be uned to define the covernge nnd renolutlon which can be achieved 
with particular orbit chnrncterlBtlcH . In Section 1 i.i below the capa- 
bilities of the fixed frequency instrument in an 800 km orbit are used ns 
a point of refe.rence for tliacusslon of the possibilities with other 
operational alternatives. This Is followed by n consideration in Section 
13.2 of how coverage will be reduced by atmospheric or surface effects. 

13.1 Global-Ocean Coverage and Resolution 

The fixed frequency MPS of Section 12.1 is designed to measure 
atmospheric pressure at the ocean surface directly below the satellite and 
on either side of the sub-satellite path at an angle of a, 14 from the 
nadir. Although the off-track scan is limited it is an Important capabi- 
lity since it allows the horizontal pressure gradient to be determined over 
the instrumental swath width. To illustrate the coverage and resolution 
obtainable with this instrument we have hosen an orbit of 796 km altitude 
at an inclination of 108° (the Seasat A orbit) . 

Figure 32 shows the surface cell size and the positions of sequential 
measurement cells near the equate" for part of an orbit. The separation 
between adjacent measurements is 200 or 240 km which is compatible with the 
1.2 mb measurement accuracy and typical pressure gradients. By comparison, 
the resolution now being demanded from temperature and water vapor sounders 
is an order of magnitude smaller. However, these variables show significant 
changes over the smaller scale while pressure gradients are generally smooth 
and between 0 and 2mb per 100 km, reaching 3 mb/ 100 km in the vicinity of 
major storms. Consequently 200 km resolution is felt to be adequate for the 
pressure sounder. 

Figure 33 shows the potential daily global coverage with an effective 
swath width of 600 km (3 measurements each separated by 200 km). The 
existing stations reporting surface pressure and typical shipping reports as 
in Figure 1 are also Indicated on this diagram. The Instrumental inte- 
gration time of 12 sec enables some 5000 measurements of surface pressure to 
be made over the oceans in 24 hours. Since these measurements are divided 
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between daytime and nighttime there are effectively 2500 new surface pressure 
data points, roughly doubling the present network. 

Current plans for an advanced operational meteorological satellite 
network call for 2 or 3 near polar arbiters. Such a multisatellite system 
could be used to Increase the density of surface pressure measurements. 

Figure 34 shows the separation of data points at various latitudes for three 
satellites with mutually displaced orbits of altitude and inclination as 
above. With these three satellites there are gaps In the coverage In the 
tropics but It Is probably adequate In view of the fairly stable atmosphere 
which prevails in this geographical belt. Coverage is more-or-less con- 
tinuous In the Important and very variable mid-latitude while the overlapping 
which occurs at high latitude could, if It were desirable, be removed and 
the resolution improved by reducing the angular scan for these parts of the 
orbit. 


It should be noted that the measurement density represented by Figure 
34 Is for a 12 hour period. The non-coincident nature of these measurements 
is characteristic of all meteorological observations from satellites and 
should not detract seriously from their potential value for numerical weather 
forecasting. 

With the standard of coverage set for one particular implementation 
of the fixed frequency instrument we can now examine briefly the Implications 
of various alternatives. Coverage and resolution with the swept frequency 
MPS In an 800 km orbit Is essentially the same. Either Instrument in the 
lower, 500 km orbit improves the cross-track resolution to 125 km but leaves 
gaps In the coverage at mld-latltudcs as well as Increasing those gaps in 
the tropics. More satellites could be used to improve the coverage and the 
resolution or tropical coverage could be selectively increased with a high 
inclination 150°) orbit. A decrease in the integration time can be used 
to Improve the along-track separation of measurement points but this must be 
at the expense of reduced accuracy. For example, a 5 sec integration makes 
the along track resolution 100 km and the r m.s. error 1.8mb with the fixed 
frequency scheme or l.Omb with the swept frequency scheme. 

Meteorological requirements will ultimately dictate which balance bet- 
ween accur. :y, resolution, coverage and cost is most appropriate for operational 
application. 



Figure "i. Separation of data points at various latitudes with a 3 satellite network. 
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13.2 Coverage Reduction by Atmospheric and Surface Effects 

The Microwave Pressure Sounder has been designed for operation over 
the open ocean. During much of the winter and spring large areas of the 
Arctic and Antarctic Oceans become covered with Ice as Illustrated by the 
seasonal variation In polar Ice cover shown In Figures 35 and 36. The Ice 
thickness grows to 2 or 3 meters so that, provided Its reflectivity la 
high enough, measurements will not be prevented by an Inadequate knowledge 
of tl reflecting surface height. Measurements by Gloersen et al (1973) 
Indicate an emisslvlty for new Ice of ''*.9 and for old ice of '>^0.8. These 
data are consistent with backscatter measurements at 13.3 GHz (AMHR, 1975, 
p. 216) and at 35 GHz (Sacklnger and Byrd, 1972) which give normalized back- 
scatter cross-sections in the range (3 to -5db) at normal incidence and 
(3 to -lOdb) at 'v» 14®. While these results suggest a reflectivity suffi- 
ciently large to enable measurements to be made In many conditions there 
Is inadequate data available on frequency variations. The technique which 
has been adopted of combining ratio measurements to account for second order 
'’'*J^i®bions has great flexibility and will probably cover this situation 
with tolerable accuracy. Occasional inadequate compliance with the design 
conditions or reduced reflectivity can be anticipated and these would lead 
to some reduction In coverage which it Is not possible to quantify at this 
time. Experimental measurements with a concept proving flight are needed to 
assure the validity of MPS measured surface pressure values In ice cover 
situations. 

The total error calculations. Tables 13 and 16, show that measurements 
will become Interrupted when cloud cover is greater than 1 km of density 
1 gm m“^ equivalent to 0.1 gm cm“2 integrated liquid water. Although global 
cloud cover Is typically near 50%, measurements by Staelln et al (1976) with 
NEMS show that the liquid water content of the atmosphere over the oceans 
exceeds .1 gm cm"^ for only a small fraction of the time. Figure 37. Some 
care must be taken in applying these results to estimate MPS coverage since 
the large angular beam width of the Instrument means that these are averages 
over 200 x 300 km ground resolution cells. Somewhat more frequent occurrence 
of such high cloud density can be expected for the 10 x 100 km cell siz.i of 
the MPS. Rain cannot be expected to reduce the coverage much further since 
It will usually be accompanied by cloud which is denser than the above limit. 
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Figure 36. Extent 


of polar Ice cover, Southern Hemisphere 
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Figure 37. Observations of liquid water content of the atmosphere 
over the oceans. Dots: one or two occurrences; cross- 

hatching: more than two (From Staelin et al, 1976). 




As noted in Section 12 some sea states will also prevent measurements 
being made but these adverse conditions occur infrequently. Also, those 
orbits for which solar radiation might also be reflected into the receiver 
aperture should also be avoided for, although this noise source is one 
which will be removed by the direct measurement, it will reduce the accuracy 
of pressure measurement. 

In contrast to these factors which reduce coverage, it may be that 
some measurements over land will be possible if the surface height is 
reasonably uniform and its reflection characteristics are appropriate. 
However, insufficient data on surface reflectivity exists to specify these 
areas and Indeed it is not a worthwhile exercise to try to do so at this 
time since the locations can more easily be identified with confidence by 
operational experience. 

It is difficult to quantify the total reduction in coverage which can 
be expected from the effects of the atmospheric and surface factors. How- 
ever it is probably not unreasonable to anticipate for the 800 km orbit a 
10% loss of data with the fixed frequency instrument and a 20% loss for the 
swept frequency instrument. For the lower, 500 km orbit these losses would 
perhaps be nearly halved. The most unfortunate aspect of this problem of 
reduced coverage is that the data would be lost from the most interesting 
meteorological features viz cyclones, with their associated obscuring heavy 
cloud and where surface pressure measurements could provide valuable infor- 
mation on the depth of the depression. Nevertheless, if technology con- 
tinues to advance at the present pace, increases in millimeter wavelength 
source powers should, within the foreseeable future, make it possible to 
measure surface pressure under these adverse conditions. 

Microwave Pressure Sounders with the design specified in Section 12 
on a three satellite network would provide some 7000 globally distributed 
surface pressure values twice daily with an accuracy of 1 to 2 mb and 
contribute significantly to the bank of meteorological data available for 
numerical weather forecasting. 
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lA. Concluslona 

Atmospheric pressure at the Earth's surface is arguably the moat 
important meteorological parameter since it is the driving force for 
atmospheric circulation and has historically served as the primary basis 
for synoptic weather forecasting. The Microwave Pressure Sounder proposed 
here would provide consistent measurements over the at present inadequately 
monitored oceans to supplement the data available from ground based stations 
and approach the requirements of the GARP program for a comprehensive 
meteorological data set. As well as contributing to the initial conditions 
for large scale numerical forecasting, the data would be of considerable 
value for atmospheric dynamic studies with General Circulation Models. 

Because a surface resolution of 100 to 200 km can be achieved, its use would 
extend to localized forecasting and to oceanography by assisting with the 
interpretation of sea state, surface wind and sea surface topography ob- 
servations from other satellite-borne sensors. Recent publication of a 
Climate Plan (ICAS, 1977) has stimulated vigorous activity in climate research. 
A significant addition to the data resources available for this program 
would be the consistent set of global pressure fields provided daily over a 
long period of time by operationally deployed Microwave Pressure Sounders. 

Absorption by a vertical atmospheric column at a frequency in the 
wing of the oxygen band centered on 60 GHz is a measure of the total amount 
of oxygen in the column. The surface pressure can be deduced from this 
measurement which must be made with an active Instrument. The selection of 
optimum operating frequencies is based upon accepted models of surface re- 
flection, oxygen, water vapor and cloud absorption. By combining the ratios 
of transmission at three pairs of frequencies the measurement of surface 
pressure is made essentially Independent of variations in the atmospheric 
temperature and humidity profiles, cloud cover and sea state. This means 
that the final output of the Instrument is a direct measurement of surface 
pressure and therefore no sophisticated inversion procedure is required. A 
numerical simulation experiment using real profiles defined by radiosonde 
observations has shotm that the error due to atmospheric profile variability 
xs less than 0.5mb. 


The set of operating frequencies used in the specific designs to 
show the potential capabilities of the MPS are not necessarily those which 
should be used in a flight experiment. Current investigations show that 
the sensitivity to pressure can be maintained or slightly Improved with 
sets which have less susceptibility to errors from broken cloud states. 

While the. effects of changes In the antenna size and Integration time are 
readily discernable from the illustrative systems, the consequences of 
changing the operating frequencies cannot In general be simply defined 
but Instead need a detailed analysis. For example, If the channel at 
52.8 GHz were changed to 53.3 GHz the other operating frequencies must be 
altered to retain temperature, water vapor and background insensitivities, 
the accuracy of the fixed frequency system would be improved by ^<20% and a 
slight loss of data would result. A better defin?*‘‘*on of accuracy, re- 
solution and coverage requirements, within the limitations presented In 
this report, will assist In determining the best set of operating frequencies. 

The choices available for hardware implementation have been Illus- 
trated by system definitions based on fixed and swept frequency alter- 
natives. Pull analyses of the errors define the balance between accuracy 
and resolution which are available with each of these designs. Coverage 
and cost are Introduced as further factors by a consideration of operational 
application alternatives. The fixed frequency Instrument in an 800 km 
near-polar orbit can provide surface pressure with an r.m.s. error of 1.2mb 
and a surface resolution of 10 x 100 km. With these measurements separated 
by 'V 200 km both across and along the sub-satellite path a three Instrument 
network would provide satisfactory global-ocean coverage at 12 hour inter- 
vals. Resolution and coverage with a swept frequency instrument would be 
similar but its accuracy Is somewhat better at ±0.6rob. This Improvement In 
performance Is however at the expense of Increased susceptibility to data 
loss through adverse sea state or cloud conditions and Increased component 
development costs. Meteorological or other user requirements will ulti- 
mately dictate the most appropriate system design and operational application 
mode. 


While some of the proposed techniques are relatively untested, the 
designs have been based, with the exception of the transmitter subsystem, 
on existing technology. The state-of-the-art for power output from railli- 
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meter wave aources is anticipated to be 2W within the next year and this will 
satisfy the requirements of the fixed frequency design. In contrast, the 
awept-frequency system will require rather more development of suitable 
sources. Sources more powerful than 2W would extend the range of atmospheric 
conditions in which the pressure sounder would operate with sufficient 
accuracy to be useful. 

Although It is recognized that these active multlfrcquency tneasure- 
tnenta Include information on parameters other than surface pressure, no 
consideration has as yet been given to the extraction of data on water 

vapor and liquid water content of the atmosphere, sea state or surface level 
which are potentially available. 

A continuing program of experiment definition studies. Instrumentation 
evaluation and ground based atmospheric transmission measurements is being 
pursued at the Jet Propulsion Laboratory in cooperation with Herlot-Watt 
University and the Appleton Laboratory in the U.K. This work is directed 
toward a proof-of-concept Shuttle/Spacelab flight in the early 1980’s in 
preparation for operational application by the late 1980's. 
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LTST OF SYMBO LS 


partial list identifies the most important symbols and those which 
used in more than one section. The section or subsection in which 
iCirst used is also given. Other symbols have been introduced and 
:n the text but their use is generally confined to one section only, 
^e convention the same symbol has occasionally been used with two 
meanings but this should not cause confusion since the correct 
31 clear from the context. 



3.2 

Antenna area 

3l,..b,_c 

3.3 

Coefficients of polynomial 

Z)_ 

11.1 

Duty Cycle 

frlf 

3.2 

Bandwidth 


3.2 

Receiver noise figure 


7.4 

Fractional accuracy 


3.2 

Height of satellite 


3.2 

Noise power 


7.4 

Noise power in channel i 

S’t.Pj.Pb 

3.2 

Transmitted, received, detected power 


4.2 

Pressure 

S?8 

3.3 

Surface pressure 

-XA' 

3.2 

Antenna radius 


7.1 

Index of surface pressure 


7.1 

Instrumental, atmospheric, background 
factors of S 


3.2 

Temperature 

1 

1 

3.2 

Integration time 


11.1 

Satellite velocity 


7.2 

Indices 

n? 

4.3 

Water vapor mass mixing ratio 


4.1 

Absorption coefficient (db km~^) 

—a p p 

A 

3.3 

Integrated absorption coefficient for 


OXVGGD. wafGr \7r»r> nr /4nf1 f i Imi- r- r'li-i a 1 


Transmit, receive antenna feed efficiency 

3db beam width 

Wavelength 

Frequency 

Operating frequency 

Normalized backscatter cross-section at 
nadir 

Atmospheric transmissivity 
Transmissivity of the atmosphere due to 
oxygen, water vapor, and other 
constituents 


